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ABSTRACT 

This  paper  documents  a  computer  program  to  be  used  in  solving 
nonlinear  programming  problems  with  nonlinear  programming  problems 
in  the  constraints.  The  program,  named  INSUMT,  is  used  with  the 
standard  program,  named  SUMT,  which  implements  the  sequential  uncon¬ 
strained  minimization  technique  for  nonlinear  programming.  SUMT 
calls  INSUMT  when  it  is  necessary  to  solve  a  nonlinear  program  in 
a  constraint.  The  INSUMT  program,  together  with  a  fairly  complete 
example  of  its  use,  is  included  in  the  documentation. 

Theory  and  applications  of  the  models  which  can  be  solved  using 
this  program  are  treated  in  two  companion  papers. 
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INTRODUCTION 


Reference  [2]  formulates  a  class  of  mathematical  programs  with 
optimization  problems  in  the  constraints.  Reference  [3]  gives 
models  of  a  number  of  defense  problems  which  can  be  approached 
using  the  techniques.  The  present  paper  documents  a  computer 
program  for  solving  mathematical  programs  with  nonlinear  program¬ 
ming  problems  in  the  constraints . 

Fiacco  and  McCormick  in  Reference  [4]  present  theory  and 
computational  aspects  of  the  sequential  unconstrained  minimiza¬ 
tion  technique  (SUMT)  for  nonlinear  programming.  Applications  of 
nonlinear  programming  to  a  number  of  practical  problems  are  presented 
in  Reference  [1].  The  initial  computer  program  for  SUMT  is  docu¬ 
mented  in  Reference  [5],  The  most  recent  computer  program  for  SUMT 
is  documented  in  Reference  [6].  It  handles  equality  constraints, 
and  includes  routines  for  computation  of  first  and  second  partial 
derivatives  using  function  values.  Various  unconstrained  minimiza¬ 
tion  methods  are  available  in  the  SUMT  program,  requiring  function 
values  and  first  partial  derivatives  or  requiring  function  values 
and  first  and  second  partial  derivatives.  Numerical  differentiation 
routines  facilitate  the  use  of  various  unconstrained  minimization 
methods.  These  are  necessary  in  solving  mathematical  programs  with 
nonlinear  programs  in  the  constraints,  since  the  partial  derivatives 
of  the  constraints  cannot  be  stated  explicitly. 

The  present  paper  is  designed  to  be  a  supplement  to  Reference 
[6]  and  to  be  used  with  it. 
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PROBLEM  DESCRIPTION  AND  SUMMARY  OF  COMPUTATIONAL  PROCEDURE 

One  of  the  mathematical  programs  considered  in  Reference  [2]  is 

to  choose  vectors  x  =  (x, ,  ....  x  )  and  v1  =  (v^,  ....  v.1  )  for 

?  rr  1?  5  k. 

i  =  1,  . . . ,  m  to 

minimize  f(x) 
x  e  X 


subject  to 

h^x)  =  min  g1(x,  v1)  >0  ,  i  =  1,  2,  . ..,  m  . 

i  „i 
v  e  V 

It  is  shown  there  that  if  g1(x,  v1)  is  concave  in  x  on  a  convex  set 
X  for  every  v1  e  V1,  then  h^(x)  is  concave  on  X  where  the  convex 
set  X  may  be  defined  by  inequality  and/or  equality  constraints.  If, 
in  addition,  f(x)  is  a  convex  function  on  X  ,  then  the  mathematical 
program  is  convex. 

To  outline  the  computational  technique,  it  is  useful  to 
differentiate  between  the  ’’outside  program,” 

min  f(x) 
x  e  X 


subject  to 


tuCx)  >  0  ,  i  =  1,  2, 


and  the  ”ith  inside  program,” 


min  g1(x,  v1) 

i  „i 
V  s  V 


m  , 
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When  convenient  the  distinction  among  the  m  inside  programs  will  be 
dropped,  considering  only  the  generic  problem: 

min  g(x,  v) 
v  e  V 

The  constraint  functions  fn(x)  in  the  outside  problem  are 
implicit  in  that  their  values  depend  upon  the  solution  of  the  inside 
problem,  which  in  turn  depends  upon  the  value  of  x  .  Thus,  a  solu¬ 
tion  technique  for  the  overall  problem  must  not  rely  on  an  explicit 
functional  form  for  fn(x). 

The  computer  program  described  in  this  paper,  called  INSUMT, 
is  based  on  SUMT.  The  standard  SUMT  program  is  used  for  the  out¬ 
side  problem.  The  new  INSUMT  program  is  used  to  solve  the  inside 
problem.  SUMT  and  INSUMT  are  iterative  routines.  Let  x  denote 
the  value  of  x  for  iteration  k  of  the  outside  problem  and 

ft 

v  denote  iteration  £  for  the  inside  problem.  The  solution 

procedure  is  initialized  by  the  user  supplying  x°  and  v°. 

During  iteration  k  for  the  outside  problem,  the  inside  problem 

is  to  choose  v  e  V  to  minimize  g(x  ,  v).  At  iteration  £  of  the 

a 

inside  problem,  there  is  a  value  v  .  This  value  is  then  used  by 

£+1 

INSUMT  to  generate  a  new  value  v  .  Continuing  in  this  manner,  the 

k 

inside  problem  is  solved,  yielding  v*(x  )  such  that 

g(xk,  v*(xk))  =  min  g(xk,  v) 
v  e  V 

Control  then  passes  to  the  outside  problem  which  generates  a  new 
k+1 

solution  value  x  .  The  process  is  repeated  until  the  sequence 
{  xk  \  converges . 

The  routine  INSUMT  resides  in  core  with  SUMT  and  may  be  called 
to  evaluate  more  than  one  constraint  of  the  outside  problem.  For 
each  such  constraint,  user-supplied  subroutines  of  INSUMT  provide 
information  concerning  the  functional  form  of  g1(x,v1)  and  of  the 
constraint  set  V1.  User-supplied  subroutines  for  SUMT  give  analogous 
information  for  f(x)  and  X  . 
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The  SUMT  and  INSUMT  programs  iteratively  generate  values  of  the 
solution  variable  by  optimizing  an  unconstrained  penalty  function. 
Specifically,  iteration  i  +  1  generates  v  from  v  and  from  the 
functional  forms  of  the  objective  function  of  the  inside  program  and 
the  constraint  set.  The  computation  calls  for  the  first,  and  some¬ 
times  second,  partial  derivatives  of  the  relevant  functions.  These 
may  be  specified  by  the  user.  However,  for  the  outside  problem, 
explicit  derivatives  of  fu(x)  cannot  be  given.  SUMT  uses  numerical 
differentiation  in  this  case  to  approximate  the  derivatives  to  be 
used  in  generating  from  x  . 
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DESCRIPTION  OF  COMPUTER  PROGRAM 


The  program  is  modular  in  structure.  It  consists  basically  of  the 
SUMT  program  (subroutines  MAIN  through  XMOVE),  one  set  of  user-supplied 
subroutines  containing  function  evaluations  for  the  outside  mathema¬ 
tical  program  (READIN,  RESTNT,  GRAD1  and  MATRIX),  the  INSUMT  program 
(subroutines  SUB  through  XMOVES),  and  one  set  of  user-supplied  sub¬ 
routines  containing  function  evaluations  for  the  one  or  more  inside 
mathematical  programs  (READIX,  RESTNX,  GRAD1S  and  MATRXX). 

The  workings  of  the  SUMT  subroutines  MAIN  through  XMOVE  are 
described  in  Reference  [6].  The  user-supplied  subroutines  for  the 
outside  program  perform  the  following  tasks  when  used  with  INSUMT. 

READIN  reads  in  the  parameters  used  in  evaluating  the  objective 
function  and  constraints. 

RESTNT  evaluates  the  objective  function  and  constraints  of  the 
outside  program  for  a  value  of  x  .  Except  where  constraints  con¬ 
tain  inside  mathematical  programs,  the  functions  are  evaluated  for 
the  current  value  of  x  .  Otherwise,  RESTNT  calls  SUB,  identifying 
the  inside  mathematical  program  to  be  solved  for  the  current  value 
of  x  ,  and  SUB  returns  the  value  v*(x).  RESTNT  then  evaluates 
h(x)  =  g(x,  v*(x)). 

GRAD1  evaluates  the  first  partial  derivatives  of  the  objective 
function  and  constraints  in  two  ways .  Where  the  constraints  do  not 
contain  inside  mathematical  programs,  for  the  current  x  the  first 
partial  derivatives  are  evaluated  either  explicitly  or  using  numeri¬ 
cal  differentiation.  Otherwise,  GRAD1  uses  DIFF1  to  numerically 
differentiate  the  constraints  which  include  inside  mathematical 
programs . 

MATRIX  evaluates  the  second  partial  derivatives  of  the  objective 
function  and  constraints  of  the  outside  mathematical  program  in  two 
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ways .  Where  the  constraints  do  not  contain  inside  mathematical 
programs,  for  the  current  x  the  second  partial  derivatives  are 
evaluated  either  explicitly  or  using  numerical  differentiation. 
Otherwise,  MATRIX  uses  DIFF2  to  numerically  differentiate  the 
constraints  which  include  inside  mathematical  programs. 

The  INSUMT  program  consists  of  subroutines  SUB  through  XMOVES. 
Three  basic  changes  are  made  in  SUMT  to  result  in  INSUMT. 

First,  all  subroutine  names  are  changed  by  simply  adding  an 
S  to  the  end  except  where  this  results  in  more  than  six  letters, 
in  which  case  the  last  two  letters  are  replaced  by  X  (e.g.,  OPT 
becomes  OPTS  and  RESTNT  becomes  RESTNX).  All  calls  to  subroutines 
within  INSUMT  are  changed  to  include  the  revised  names  of  the 
called  subroutines.  SUMT  calls  INSUMT  only  by  RESTNT  calling  SUB. 
After  SUB  is  called  INSUMT  calls  only  the  subroutines  within  INSUMT 
until  SUB  returns  to  RESTNT. 

Second,  all  labeled  COMMON  arrays  are  changed  by  adding  an  S 
to  the  end,  except  where  this  results  in  more  than  six  letters,  in 
which  case  the  last  two  letters  are  replaced  by  X  (e.g.,  SHARE 
becomes  SHARES  and  CONPAR  becomes  CONPAX). 

Third,  subroutine  MAIN  is  modified  to  obtain  subroutine  SUB  by 
changing  the  name  from  MAIN  to  SUB,  and  by  modifying  the  subroutine 
to  read  in  data  for  the  inside  mathematical  programs  the  first  time 
each  problem  is  solved  and  to  save  the  data  in  an  array  for  subse¬ 
quent  solutions  of  the  problems. 

User-supplied  subroutines  READIX,  RESTNX,  GRAD1S  and  MATRXX 
are  called  by  INSUMT.  Depending  on  a  parameter  denoting  the  inside 
mathematical  program  being  solved,  the  subroutines  calculate  appro¬ 
priate  function  values,  first  partial  derivatives,  or  second  partial 
derivatives . 

In  the  first  several  runs  of  a  program,  it  is  necessary  to 
observe  the  intermediate  points  of  the  solution  of  the  inside 
programs.  RESTNT  thus  prints  out  which  inside  program  is  being 
solved,  and  INSUMT  prints  out  the  intermediate  points.  After  it 
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is  established  that  the  inside  programs  are  being  successfully  solved, 
printing  of  their  points  is  suppressed  by  deleting  the  print  statements 
from  RESTNT  and  by  modifying  BODYS,  CONVRX,  ESTIMS,  FEASS,  INVERX, 

OPTS,  OUTPUX,  PUNCHS,  TECHECX  and  TIMECS .  The  modifications  are 
described  later. 
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EXAMPLE  PROBLEM 


An  example  problem  is  used  to  describe  the  use  of  INSUMT  with 
SUMT. 

The  problem  is  to  choose  x^,  xn  to 

minimize  xn  +  ...+  x 
1  n 


subject  to 


minimum  xi^vi 

vn  +. . .+  v  <  n 

1  n  — 


2>  +...+  vv2)2 


2  rl 


5  2  •  5  2 

minimum  x*  (v  -  2)  +...+  xn  (vn  -  2) 

vf  + . . .+  vz  <  n 
1  n  — 


£  r2 


To  illustrate  the  problem,  take  n  =  4,  r^  =  4,  r^  =  4  and  let  the 
starting  point  be  x1  =  x2  =  x^  =  x4  =  4  for  the  outside  program  and 

=  .5  for  both  inside  programs .  The  value  of  the 
objective  function  of  the  outside  program  is 

X1  +  x2  +  x3  +  x4  =  16  . 


The  objective  function  of  the  first  inside  program  is 

4(.5  -  2 )2  +  4( . 5  -  2)2  +  4( . 5  -  2)2  +  4( . 5  -  2)2  =  36  , 

and  since  .5+  .5+  .5+  .5=2<4  the  first  inside  program  starting 
point  is  feasible.  For  the  second  inside  program  the  objective 
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function  is 


2( . 5  -  2)2  +  2( . 5  -  2 )2  +  2( . 5  -  2)2  +  2(.5  -  2)2  =  18  , 

2  2  2  2 

and  since  .5  +.5  +.5  +.5  =  1  <  4  the  second  inside  program 

starting  point  is  feasible.  Since  36  >  r^  (=  4),  18  >  r^  (=  4)  , 

the  three  starting  points,  one  outside  and  two  inside,  provide  a 
feasible  point  for  all  three  programs. 

The  solution  to  the  example  is  =  1,  with 

Vf  =  =  V3  =  =  1  in  both  inside  programs,  yielding 

xl+x2+x2+x4=l+l+l+l=4 


and 


1(1  -  2)2  +  1(1  -  2)2  +  1(1  -  2)2  +  1(1  -  2)2  =  4 

s.t. 

(1  +1+1  +  1  =  4)  <4 


>  4 


9  2  2 

1(1  -  2;  +  r(l 

s.t. 

..2  n2  . 2  .2 
(1  +1  +1  +1 


2)2  +  12(  1  -  2)2  +  12(1  -  2)2 
4)  <  4 


4 


>  4 
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COMPUTER  PROGRAM  INCLUDING  USER-SUPPLIED  SUBROUTINES  FOR 
EXAMPLE  PROBLEM 

In  this  section  the  user-supplied  subroutines  for  the  outside 
program  are  presented,  followed  by  the  INSUMT  program,  followed 
by  the  user-supplied  subroutines  for  the  inside  programs.  The 
SUMT  program  is  not  supplied,  being  documented  in  Reference  [6]. 

RESTNT  contains  print  statements  for  each  call  for  solution 
and  end  of  solution  of  an  inside  program.  The  INSUMT  subroutines 
SUB,  BODYS,  CHCXEX,  CONVRX,  ESTIMS,  FEASS,  INVERX,  OPTS,  OUTPUX, 
PUNCHS,  TCHECX  and  TIMECS  contain  print  statements  for  the  points 
of  the  inside  programs.  Thus  the  program  as  listed  prints  all 
intermediate  points.  However,  the  changes  necessary  to  suppress 
printing  are  indicated  in  heavy  markings  directly  on  the  listing. 
Simply  remove  the  boxed  statements,  and  replace  them  by  the  state¬ 
ments  written,  if  any.  Printing  is  done  by  SUB  only  the  first  time 
it  is  called,  so  SUB  is  not  modified.  Printing  is  done  by  CHCKEX 
only  if  control  cards  dictate,  so  it  is  not  modified.  All  other 
print  statements  are  removed  to  suppress  printing  of  inside  programs. 

It  should  be  noted  that  all  data  of  the  example  problem  are 
contained  in  the  set  of  user-supplied  subroutines  for  the  outside 
program  (READIN,  RESTNT,  GRAD1,  and  MATRIX)  and  the  set  of  user- 
supplied  subroutines  for  the  inside  programs  (READIX,  RESTNX,  GRAD1S 
and  MATRXX).  For  most  problems  READIN  should  be  used  to  read  data 
for  outside  and  inside  programs,  since  READIX  is  called  each  time  an 
inside  program  is  solved,  and  it  should  not  be  used. 

Control  cards  are  read  in  the  following  order:  Outside  program, 
first  inside  program,  second  inside  program,  and  so  on  if  there  are 
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more  inside  programs.  SUB  saves  the  control  card  data,  and  inside 
program  starting  points,  and  provides  these  data  each  time  an 
inside  program  is  solved.  Dimensions  on  PARS  are  presently  (2,47), 
including  up  to  20  variables  in  the  inside  program.  This  would 
have  to  be  changed  for  more  than  2  inside  programs  or  20  variables. 
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User-Supplied  Subroutines  for  Outside  Program 
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000002 

000002 

000003 


SUBROUTINE  REAfllN 
9999  CONTINUE 
RETURN 
END 


000005 

000005 

000005 

000005 

000006 

000007 

000010 

00001? 

000016 

000016 

88881! 

000023 

OOOO?4 

000027 

000031 

000036 

000041 

000041 

000042 

8888:: 

000051 

000057 

000061 

8888^ 

000064 

0000^5 

000070 

OOOO72 

000076 

000 1 02 
000102 

000103 

000107 

000107 

88811$ 

000u6 

000126 

000l?6 

000127 


SUBROUTINE  RESTNT ( IN.VAL) 


i* 


COMMON/PROQ/ISP 
COMMON/OUT/NOUT»XT (100)fV(100) 

COMMON/$HARt/  X(100)t  OEL<l00)f  A ( 1 0 0 ♦ 1 0 0 ) • N. M*  MN*NPl*NMl 
fn.n 

VAL*0* 

IP  (IN) 100*l00ip00 
DO  150  j*  1  *  N 
val  .  vAl  *  X ( j) 

60  TO  9999 

200  IP  <  IN— 1  )  300*300*4-00 
300  I$P*1 
NOUT*N 

00  3°5  J*1*N 

,  IF(X(j)-o.)joS3oS3o5 

30*  X(j)b1,e-Ob 
305  xT  (jj  ■  x  (  |) 

- »»rETr,367) 


307  format  1 1 1  Hi  PROBLEM  a> 

- Call  5Ub  4 - 


WRITE (6*308) 


F0RmAT(18m1  PRflBl  EM  A  SOLVED) 
D0^7TrjiTrTN 
3?0  VAl  .  VA,  «  X(j)*(V(j)-?.)»»2 
VAL-VAL-*. 


*00 


GO  TO  9999 
I  SP»j 
NOUT  »N 


OO  *05  Jb 1  * N 
IF(X(j>-0.)40*,404*405 
404  X(j)b1.E-0b 

4°5  XT(J)»X (j) _ 

WHITE (6. *b7) 

407  FORMAT ( 1 1  Hi  PROBLEM  B) 


call  SUB 

WRITE<6*40S> 

406  PORHAT ( 1 ftHl  PROBLEM  P  SOLVED) 

DO  420  J,1#N 

420  VAL  •  VAL  ♦  X ( J) • 5  •  (V(J)-2#)«#2 
VAl*vAl-4# 

GO  TO  9999 
9999  CONTINUE 

return 

END 


10000070 


1? 


000003 

oooooj* 

0000O5 

000010 

50 

000011 

loo 

000013 

000017 

150 

nmi 

200 

000022 

000023 

9999 

surroutine  gradkin) 
common/share/  x ( 1 00) , 

00  50  J«1 ,N 

del  < *5. 

IF (IN) 100,100,200 
DO  i'll  j»l  »N 

■  1. 

GO  T0  9999 
CALL  DlFFi (IN) 

GO  TO  9999 

RETURN 

end 


OEL(100) t 


A(100#lOO,#NfMf  MN  »  NP 1 • NM 1 


10000070 


IB 


000005 

000005 

000006 

000007 

000010 

000011 

000013 

000013 

000014 


SUBROUTINE  HATPIXUN.U 

COMMON/SHARE^  x'lOO1*  OEL<100>»  A  ‘  100 *  100 > *N.M*  MNtNPltNMl 
IF (IN) lOOtlOOtZOO 
lOO  L*1 

G°  T°  9999 
200  CALL  DIFfZ'IN* 

GO  TO  9999 
9999  CONTINUE 

return 

end 


'»sf^TtwjrjTr 


10000070 
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INSUMT 
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subroutine  sub 


000002 

000002 

ooooo2 

000002 

000002 

000002 

000002 

000002 


000002 

000002 

000002 

000002 

nml 


0000 1 2 
000036 

888815 

OOOO54 

000055 

00005ft 

OOOloo 

000112 

000136 

000142 

000172 

000202 

000206 

000216 

OOO??6 

000232 

000242 

000244 

000?46 

000247 

000251 


C 

C  MARCH  1971 

C 

C  MAIN  IS  THE  PROGRAM  THAT  INITIATES  THE  SUMT  ALGORITHM*  THE  INPUT  OF 
C  PARAMETERS*  OPTIONS*  AND  STARTING  POINT  IS  DONE  IN  MAIN.  AFTER  ThE 
C  SOLUTION  OF  ONE  NLP  PROBLEM  MaIN  LOOKS  FO®  OAtA  port  ANOTHER  NLP  PRO®. 

c 

COMMON/PROR/ISP 
COMMON/OUT/NOUT . XT ( I 00>  * v  U 00) 

COMM0N/In/W(100> 


c 


COMMON/ShARES/X (100) ,  DEL  ( 1 00 ) *  A < 1 00* I  00 > *N , M ,  MN * NP I  * NM I 
COMMON  /EQALS/H*  H  j  f  MZ 

COMMON  /OPTNSS/NTl *NT2*NT3fNT4*NT5*NT6*NT7tNTfl*NT9,NTlO 
COMMON  /vaLUES/F,G*^0*rSIgmA#  Rj ( 200 ) *  Rh° 

COMMON/CRSTS/OELx (1001 *  DELX0(100)v  RM01N*RATI0»  EPSI*  THETAo* 

1  HSlGif  GI*  Xi  (100)*  X2doo)*  X**dOO)*  XRpdOO)*  XR]dOO)t  PRl  t 

2  PR2 * P ) *  pA*  RjI(200)*  DOTT*  PGRAp(IoO)*  DlAftPOO)* 

3  RrEV3,A0ELx#  ^TcTR,  NU*INI,  ^RHASe*  nSATI$ 

COMMON/EXPOPX  /  NEX0P1*  Np X0P2 *  XEPI*  XEP2 


C 


C 


C 

C 

10 

C 


0  ImENS I  ON  1FTS(2) ,PARS(2.47) 
DATA  IFTS/o*0/ 


5 


00  5  J-l.NOUT 


W (j)*XT ( j) 
IF ( 1 F  TS ( 1 SP 


,,10*10,20 


PARAMETER  CARO 

REAO  (5*50>  EPSI *RHOIN,THETAO*RATIO*TMMAX*M,NtMZ 

INITIAL  x  VERCToR  card  format 
REAO  (5,60,  (X( I, , 1*1 ,N, 

NTCTR*0 
NP 1 *N  ♦ * 
nmUn-1 


c  SUBROUTINE  ReADIN  is  UNDER  PROGRAMMER  CONTROL 
CALL  REA01X 

C  OPTION  CARD  FOLLOWS  PROGRAMMERS  OaTa 

REAO  (5*8?)  NTi*NT2*NT3*NT4»NT5»NT6*NT7*NTb*NT9»NII0 
WR 1 T  p  (6,aA0) 

WRITE  (6,1201  N,M,MZ*TMMAX*RHOlNfRATiO,EPSi, THETAO 
WRITE  6*130 

WRITE  (6,80)  NT1 *NT2*NT3*NT4*NT5*NT6*NT7*NT8*NT9*NTI0 
C— READ  TOLERANCES 

read  (5,6o,  xEpi,  xep2 
WRITE  <6,90> 


WRITE  (6*70)  XpP],  XEP2 
C— READ  SECOND  OPTION  CARD 

REAO  (5*8°) *NEXOPl *NEX0P2 
WRITE  (6,100) 

wRiTe  (6,80)  NpxOPj ,NExOP2 
PARS  (ISP*  P  ■  EPSI 
PARS (ISP,  2)  *  RhO I N 
PARS  (ISP,  ->  *  THETAO 

pars  (iSP,  4,  «  Ratio 

PARS  (ISP*  5)  ■  TMMAX 


000030 

000040 

000050 

000060 

000070 

0000*0 


000 1  BO 

000190 

00025q 

000260 

88812! 

0002q0 

OOO3OO 

000320 

OOO33O 

0003*0 

000350 

000360 

000370 

0003Bo 

000390 

000400 

000410 

0004^0 

000430 

000440 

000450 


22 


000252 

00025* 

000256 

000260 

000262 

00026* 

000271 

000273 

000275 

000277 

00030l 

000303 

00030* 

000307 

000311 

000313 

000315 

000317 

00032o 

000322 

000323 

000325 

000327 

000331 

000333 

ooo335 

000337 

000337 

gggSti 

0003** 

0003*6 

0003*7 

000351 

000353 

000355 

000357 

mm 

00037  o 
000372 
00037* 
000376 
000400 
000*02 
000*0* 
000*06 
000*10 
000*12 
000*1* 
000*15 
000*17 
000*20 
000*22 

!SS*St 

000*30 


PARS  (ISP.  6) 

■  M 

PARS  (ISP.  7) 

■  N 

PARS  (ISP,  8) 

■  M2 

DO  *  K»1,20 

L«K*r 

PARS ( I SP ,L)  ■ 

X  ( K ) 

P*RS  (ISp.29) 
P*RS  ( I Sp * 3y » 
P*RS  ( I Sp  *3I) 
P*RS  ( I Spi 32) 
P*RS  <lSp*33> 
P*RS  ( I Sp» 3*) 
P*RS  (ISP|35) 


P*R; 

PARS 

PARS 

PARS 


(ISP.36> 

(PP.37 


<ISP*38 
(  ISP, 39) 
(lSP,*0> 
PARS  (ISP.*I) 
PARS  <ISP.*2) 
PARS  ( i SP» *3 » 
(ISP.**) 
(ISP,*5> 
'ISP'46’ 
(ISP. *7* 


PARS 

Paps 

PARS 

pars 


iets<isp).i 

G°  TO  25 


NTl 

nT2 

NT* 

NT* 

nTI 

nt6 

nTr 

NT10 
■  NTn 
X£PI 
XEP2 
XEP3 
NeXOPI 

BII8S1 

NEX0P4 

NEXOP5 


PAR$<I$PiI> 

PaR$<ISP*2> 

■  pARS<l$Ptf) 

■  PARS ( I Sp  *  ^) 

*  pars jj sp • *) 
P*RS(1SP»”) 

P*RS( ISP.7) 


20  EPSI  ■ 

RHOIN« 

TheTAO 
B*TI0 
ThmAx 
M 

N  -  - - -  -  . 

mZ«  P*RS(ISP»8) 

0°  21  Kal ,20 
L"  K*8 

21  XIK)  ■  PARS ( I SP ’L) 
NT1  ■  P*RS ( ISP, 29) 

PARS ( ISP. 3o) 
PARS ( lSP.31 ) 
Par5 1 1 Sp,32  > 
PARS<ISP*33> 
PABSdSP.34) 
PARSdSP.35) 
PARSIlSP.36) 
PARS ( I SP . 3  7 ) 
PARS! ISP.38) 
PARS ( ISP.39) 
PARS ( iSP.*0) 
PARSIlSP, *1) 
_  PARSIlSP.*2) 
NEXopl.pA«SMSPt43) 

NExOP2“PARS'ISP*44» 
NEX0P3"P»R$  < ISP. 45) 
nEXOP*«P*RS(ISP.*6) 
NEX0P5«P*RS(ISP»*7> 


NT2 

NT3 

N74 

NT5 

NT  6 

NT7 

NT8 

nT9 

NT  1 0 

NT1I 

XEPl 

XEP2 

XEP3 


000432 

60  TO  25 

000433 

c 

25 

CALL  SETs ( TMMAX ) 

000435 

CALL  TIMECS 

000436 

nphase** 

888*e8 

00^437 

C - 

JUST  TO  GET  AN  INITIAL  printout 

CALL  EVALUS 

000440 

P0=0.0 

000500 

OOO44I 

G*0, 0 

OOO5IO 

000442 

HsO.O 

OOO52O 

000443 

RSIGMA*0, 0 

OOO53O 

000444 

CALL  OUTPUX  (2) 

000445 

CALL  stores 

000446 

00045“ 

IF  (NEXOPl ,GT , 1 )  CALL  CHCKEX 

IF  (NEXOPI.EQ.:’)  STOP  01072 

000570 

000456 

IF  (NEXOpI »EQ«5)  STOP  01,04 

0005fl0 

000462 

call  feass 

000463 

c  NPHASE»5  IS  USEO  to  INDICATE  NO  feasible  point  exist 

GO  To  (30, 3°, 30,30, 10) ,  nphaSf 

888|?8 

000474 

30 

NPhaSE*2 

^TcTRrO 

00062o 

000475 

000630 

000476 

call  bodys 

000477 

35 

00  35  J* I • N 

000^01 

V(j)»X(j) 

000505 

c 

retuRn 

000670 

c 

PARAMETER  CARO 

000680 

000505 

50 

FORMAT  (5EU.0OU) 

00069q 

c 

INITIAL  X  VERcTOH  CARO  FORMAT 

000700 

000505 

OOO5O5 

70 

FORMAT  ( 6E 1 2 • 6) 

FORMAT  (6f20#7) 

000710 

OOO72O 

c 

OPTION  CARD  FORMAT 

OOO73O 

000505 

eo 

FORMAT  (1017) 

000750 

000505 

oO 

format {*3h0  tolerances  ) 

OOO76O 

OOO5O5 

loo 

FORMAT  (26HJ  SECONO  SET  Of  OPTIONS  ) 

FORMAT  (56h1  NONLINEAR  PROGRAMMING  RoUTiNE-SUMT  VERSION  4  SUB 

OOO77O 

OOO5O5 

1 10 

OOO70O 

ooogOo 

OOO5O5 

120 

1  ) 

FORMAT  ( IH°«5Xf 2HN*I3t6xt2HM«I3«6x*3HMZ*I3//BX« l0HMAX.  TIME*EU*7* 

I4X,2HR*EI4. - t4X*6HRATlOaEI4.7f6Xf 0hEPSILOn*E14.7*4X*6hTHETa*E14.7) 
F°RMaT  (  \ RHO  OPT i ONS  SfLECtED) 

00081 0 

000505 

130 

000820 

000505 

end 

000830 

24 


000002 

000002 

00000? 

00000? 


000002 

ooooo2 

000003 

000004 

000005 

000006 

000007 

000010 

000011 

00001? 

000020 

00002? 

0000?3 


000024 


000045 

000053 

000054 

OOOO56 

000060 

000061 

000070 

000071 


000100 
00010? 
000110 
000 1 1 6 
000124 
000125 


000127 

000130 

000133 

000142 


SUBROUTINE  BODyS 

OCTOBER  1970 

BODY  COOROINATES  THE  FLOW  AMONG  THE  SUBROUTINES  THAT  ACTUALLY  00  THE 
CALCULATIONS  REGUIREO  8Y  The  VARIOUS  PARTS  of  the  ALGORITHM. 
C0MM0N/ShARLS/X (100) f  0EUl00>*  A(l00*l°0) iNtM*  MN*NPl*NMl 
COMMON  /0PTfcSS/NTlfNT2tNT3*NT4iNT5*NT6*NT7tNTRfNT9tNTl0 
COMMON  /VALUES/F^GtPOfRSlGMA,  R J (200 ) *  rhO 

common'crsts'oelx <100> *  OELXO(lOO)*  rmoin*ratio*  EPSI*  thetao* 

1  RSlGlt  Git  Xl(100)«  X2(100>t  X3  ( 10  0 ) •  XR2(100>t  XRl(lOO).  PRl . 

2  PR2 .Pi.  F1*  Rj1(2oO)*  OOTTt  PGRAO(lOO)*  OlAG(lO0>* 

3  PREV3fAoELx*  NTcTR.  NUm INI t  NPhASE*  NSATjS 
COMM£N/CONP*X  /  NF 1 t  NF2tNF3 

nF3.2 

MN*0 

NUMINI-0 

C  OPTION  Of  GEt7iNg  INiTIaL  RhO 

CALL  RHOCOX 
CALL  EVALUS 
10  CALL  XMOVES 

GO  TO  (30*20) •  NT3 
20  CAi l  IIMFCS 

CALL  WUX  (1) 

GO  TO  40 
30  CALL  TCHECX 

C  IN  FEASIBILITY  PHASE  4  MEANS  FEAS  ACHIEVEO 
40  GO  TO  (50  .5**50*200) «  NSATIS 

50  CALl  CONVRX  (Nl) 

GO  TO  (60,10,130) ,  N1 
C  MINIMUM  ACHIEVEO  IF  Nl-1 
60 
70 


C  — 
8° 


90 

C 

C 

c 

100 


11° 
1*0 
130 
140 
C 
C 


GO  TO  (7o.80) t  nT3 
CALL  TIMECS 
call  OUTPUX  (1) 

NUMBER  Of  MINIMA  ACHIEVEO  INCREASED  by  1 

NUMINI-NUMINI^I 
MN«® 

GO  TO  (190t90t9°).  NPHASE 
CALL  ESTiMS 

finaL  Wight  have  been  caLLeo  by  eStih„c°nveRgeo  if  n2.i 
GO  TO  (l00*110*120>  *  NTA 
NT*»i  FINAL  CONVERGENCE  ON  o  ORDER  ESTlH 
ORDER  ESTIMATES.  NT*-3  converge  on  SECONO 
CALL  FINaLS(Nfi) 

*130*1*0*  ’ 

(130.140). 

( 1 30 . 1 4 


'0) 


NF  1 
NFJ 
NF- 


000850 

000860 

000870 

000880 

OOOgQO 


ESTIMATES.  NTa»2  CONVERGE  on  firs 
OROER  ESTIMATES. 


GO  TO 
GO  TO 

SeTuSn' 

RHO"RHO/RATIO 

USING  PREVIOSLY  COmPUTEO  VALUES  FOR  F.  ANO  RJ  P  IS  RECOMPUTED  WITH  THE 
NEW  V*LUE  OF  RHO. 

call  pevalx  ,  , 

C  A  VECTOR  IS  left  in  DELX(i>  by  ESTIM 
IF  (NUMINI»£)  10,150,150 
GO  TO  (10.1?0,160> .  NT7 
HAnS, 


1*0 

160 


000980 

000990 

001000 

881818 


001060 


001090 

001110 

001120 

001 1  AO 
0011*0 
001160 


OOU9O 

001200 

001210 

001220 

001240 

001250 

umi 

001290 

001300 

001310 

001320 

001330 

001340 

00l350 

001370 

001380 

OOl39o 


CALL  G"Ao5 (2) 


25 


000144 

000145 

CALL  OPTS 

GO  TO  1  lflOilTfl) ,«  NT3 

001420 

0001 53 

0  00 1 57 

WRITE  I6,?1!)) 

CALL  OUTpUX  ( 1 > 

170  CONTINUE 

001450 

0001*1 

1 80 

GO  TO  50 

000162 

190 

IF  (G)  9o» 90*200 

001460 

0001 ft4 

200 

RETURN 

0  mEans  no  feasible  point  exists 

OOl47o 

C  --- 

C 

DUAL  VALUE  GREATER  ThaN 

0014RO 

001490 

000165 

0001*5 

210 

FORMAT  <6X,30HMOVED  ON 
ENO 

extrapolation  vector  ) 

001500 

OOI5IO 

26 


subroutine  chckex 


00000? 

ooooo? 

ooooo? 

000005 

000006 

000010 

000012 

000013 

000015 

0000?? 

88883) 

041 


88! 


00005* 

10 

c.. 

888811 

oooo65 

000067 

000071 

88882f 

oodoi 

000102 

20 

000110 

30 

000112 

OOOl” 

000116 

888111 

000123 

000126 

000127 

000137 

0001*0 

40 

000143 

50 

000145 

0001*' 

60 

000150 

000153 

70 

000156 

000163 

000164 

80 

000203 

90 

000206 

000?07 

000210 

100 

000216 

000220 

no 

MARCH  1971 

CHCKER  COMPUTES  AND  LIST  THE  FIRST  PARTIAL  DERIVATIVES  USING  GRAD1 
AND  then  USING  numerical  DIFFERENCING  (DIFFi ) .  IF  REQUESTED  THE 
SeC°Nd  paRtIaL  UeRivAtIvEs  ArE  C°MrUTed  AnD  Listed  USiN(5  MaTRiX  *Nd 
diff2, 

C0MM0N/ShaRES/X(100) #  OEL(IOO)*  A(100*100) #N#Mf  MN*NPl*NMl 
COMMON  /eQALS/h*  Hi#  MZ 

D°  5  J.l," 

0EL(J*»1. 23*5678 
CONTINUE 
DO  10  I »i t  MMZ 
It!*!”1 

WRI*E  (6,170)  IN 
CALL  GRADlS'lN) 

WRITE  (6tl80>  (DEL ( J) »  J»J »N) 

C*LL  DIFFJS(IN) 

WRiTe  (6,1b5!)  ( DEL  <  J>  •J«1»N> 


IF 


GO  TO  160 


TO  *0 


(NEXOPi .LT.4) 

WRiTF  <6. *9?) 

DO  150  I ■ 1 » MHZ 
INaI-l 

WRITE  *6* 170  IN 
ITa? 

D0  30  Kal,N 
DO  ?0  ja^N 
A(K,J)bO# 

CONTINUE 

CALL  MiTRXX ( IN, IT) 

IE  (IT.EO.1)  GO  TO  l50 
D°  50  K.2,n 
KMlaK-1 
DO  *°  ja I • KM 1 
if  AS  <K* J) .tO. 0.0)  GO 
^Ex^l-S  * 

WRITE  <6,210>  K.J 
GO  TO  60 
CONTINUE 
CONTINUE 
DO  90  K« 1 ,N 

DO  70  J»K • N  „  . 

IF  (A(K,J) .Ne.o.O)  6°  TO  80 
CONTINUE 
WRITE  (6,220)  K 
GO  TO  ’0 
WRlTE  (6,200)  K, (A(K,J> ,J»1,N) 
CONTINUE 
DO  HO  Ka  1  ,N 
DO  100  jal.N 
A (K.J) »0 
CONTINUE 

WRITE (6* 1 15*  lN 


OOI53O 

0015*0 

001550 

001560 


8811^8 

OOI59O 


00 1 63q 
001640 
001650 
001660 

001670 

001680 

001690 

001710 


OOI73O 

001**0 

001750 


001*80 

001790 

001800 


88l|18 

00 1830 
001840 
001850 


OOlBZo 

001880 


00l9l0 


001930 

0019*0 

001950 

001960 

001970 

OOl’Bo 

001990 

002000 

002010 

00?0?0 

002030 

0020*0 

002050 

OOZO&O 

002070 

002080 

002090 


27 


000226 

115 

format  (13H0  CALL  DIFF2(,i2,1h)  ) 

000226 

000230 

C*LL  D1Ff2s(1N) 

DO  **°  tN 

000^32 

00  120  J*K • N 

000233 

IF  (  A  ( K  *  J )  .NE.O)  G°  TO 

130 

000236 

120 

CONTINUE 

000241 

000246 

WN I TF  (6*220)  K 

GO  TO  140 

000247 

130 

*HITF  (6#200)  K  * {A ( K  .  J ) 

*  U* 1 • ^ ) 

000266 

140 

CONTINUE 

000271 

150 

CONTINUE 

000274 

160 

CONTINUE 

000274 

r 

RETURN 

000275 

000275 

L 

170 

180 

FORMAT  (  2AH0CHECKFR. 

....constraint  NO.  •  13) 

FORMAT  ( 1H0 t24HCHECKER • 

....  1ST  PARTlAUS/(lXtE20.8*E20,8*E20.8*E20. 

16*^2q#R#^2o*8) ) 

,2nD  PaRTIALS) 

000275 

190 

F°RmAT  ( !H0t24HCHECKFR , 

000275 

000?75 

m 

FORMAT  (4Ho«OW*  13 

/(lX*E?o.fl,P2Q»P,E20»P,E20*8*P20-8’E2o«fi)  > 

FORMAT  (3h  A<*I2*1H*«12 

•  1  Oh)  .nf. 

000275 

22  0 

FORMAT  («H  R0W*l3,llH  All  ZEROS.) 

000275 

END 

002100 

002120 
00?l3o 
002140 
002150 
0021  AO 
00?170 
002180 
002190 
002200 
002210 
002220 
002230 
002240 
002250 
00??60 
002270 

S81i”2 

002300 

002310 
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SUBROUTINE  CONVRX  (Nl) 


000003 

ooooo3 

000003 

000003 


000003 

ooooo; 

000003 
000004 
00001  1 
0000?0 
000024 
000024 
000032 


oooo43 

000047 

000053 

000056 

000060 

000061 

OOOO*5 


000067 

8888;? 

000073 

000073 

000073 

OOOO^3 


0ct08ER  1970 

after  each  iteration  of  the  algorithm  to  locate  the  minimum  of  the 
PENALTY  FUNCTION,  cONVRG  DETERMINES  if  The  current  point  is  close 
enough  to  the  point  GIVING  the  minimum  VALUE  of  the  p  function, 

ni  Set  equal  to  ]  mjNjmum  Has  been  f°uNo# 

Nl  SET  EQUAL  TO  2  IF  MINIMUM  HAS  NOT  BEEN  FOUNO  ANO  TIME  IS  NOT 

Nl  set  equal  to  3  OTHERWISE 

oOTT  Set  eGuaL  to  (del  P) (INVeRSe<DEL<PEL  P)))(OEL  p)  IN  OPT 
COMMON/ShAReS/X ( ^0) *  OEL(l00>*  A  <  A • 3 ® 0 ) tNtM*  MNtNP^NM1 
COMMON  /OPTNSS/NTl tNT2tNT3tNT4tNT5fNT6tNT7,NT8fNT9,NTlU 
CflMMQM  /VA^UeS/Fi G»PO*RSlGMAf  Rj(200)»^R|jO 


10 

20 

♦38 

50 


60 


C 

C 

70 

75 


UP 


C0mm0N/CrsTs/OeLx(100) f  OELXO  < 1 0  0 ) ♦  RhOjN, RAT lO#  ePSI,  THETAOt 
1  PSIGIJ  G 1  *  ?  •  X2U00)f  x3  { lflO>  *  XR2<lOQ>f  XRi(100>’  P«1 

\  P R2iPX»  F1  •  RjM200)*  OOTT t  PGRAo(A00>t  DIAg(I00>* 

3  PReV3*Aqelx»  NTcTRt  nUm IN  I »  nPhASe*  nSATtS 

COMMON/EXPOHX  /  NE^OPd  NExOP2f  XeP*.  XEP* 

C^M^ON  /TSWS/NSWW 

(MN.lp.l)  Ql*pO 
TO  (lu#2V,30) f  nT9 
(  A0S  (  qOTT  )  ,  j_T  •  ePS  i  )  gO  TO  7o 
T°  40 

<ABS(OOTT>  ,LT,  <P1-P0>'5.0>  GO  TO  70 
TO  40 

(AOELX.LT.EPSD  GO  TO  70 
TO  (50,6°),  NSWW 
(MN.Le.1^  rETuRn 
IF  <P0*XEP2  .LT.  Ql)  G°  TO  75 
IwBtTe  (6»'6'nTf 

GO  TO  70 

call  punchs 

I  WRITE  (6.90H 

N1*-5 


IF 

GO 

S5 

IF 

GO 

IF 

GO 

IF 


FOUNO  THE 
RETURN 

Q1  -  P° 

Return 


MINIMUM  TO  THE  SUBPROBLEM* 


c 

80  FORMAT  (lOOH  APPARENTLY  ROUNOOFF  ERRORS  PREVENT  A  MORE  ACCURATE  DE 

iterminaTion  of  The  minimum  of  This  subproblem.) 

9°  FORMAT  (48h0*#**  time  IS  UP*  CALLING  EXIT  FROM  CONVRG,  ••••) 

ENO 


002330 

002340 

002350 

002360 

00?370 

0023fiQ 

002390 

002400 

002410 

002420 


002520 

002530 

002540 

00255q 

002560 

002570 

002580 

881188 

00261q 

002620 

002640 


00^67  o 
002680 
002690 
002700 

88i;i8 

003730 


003740 

002750 

003760 

002770 


00378o 
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SURROUTlNE  OIFFIS(IN) 


C 

C 

FEBIJARY  1971 

oozeoo 
00281 0 

C 

c  1 

DlFFl  COMPUTES  THE  FIRST  DERIVATIVES  RY  NUMERICAL  DIFFERENCING. 

002820 

002830 

c 

0 

-USER  CAN  CALL  FOR  DIFFERENCING  OF  SELFCTEO  FUNCTIONS 

002g*0 

002850 

000003 

000003 

000003 

000003 

COMMON/ShaReS/X(100> *  DEL ( 100) .  *(100*100) ,N.M’  MN«NPl*NMl 

CQMMQN/EXPOPX  /  NEXOP1*  NExOP2*  XEPl.  XEP2 

cOmmON/  STiRxS/xSTR,  100  )  ,  xSSS(  1 00 ) •  ODLL(IOO) 

DO  lo  J* 1  * N 

002890 

ooooo5 

10 

XSTR(j)*X(j) 

D°  30  J*l,n 

002900 

000011 

002910 

000012 

IF  (J.EO.l)  GO  TO  20 

002920 

000014 

JMl*j-l 

002930 

000015 

X(jM»)*XSTR(JM1) 

OO294O 

OOOOI7 

20 

X  (j)  *XSTr  (j)  *X|rpl 

002950 

000022 

OOOO24 

CALL  ReSTNX  ( in* ZZ2) 

X(j)*XSTR(j)-Xcpl 

CALL  RESTnx  ( I”* ZZl ) 

DEL<J,*<ZZ2-ZZi >'<2.*XEPll 

00297o 

000027 

000032 

30 

002990 

000042 

X(N)*XSTR(N) 

003000 

000043 

return 

003010 

000044 

end 

003020 

30 


SUBROUTINE  D1EE2SI1N) 


C 

003040 

C 

OCTOBER  1970 

003050 

C 

003060 

C  1 

■)IFF2  computes  the  SECONn  DERIVATIVES  by  NUMERICAL  DIFFERENCING. 

003070 

C 

0030g0 

000003 

COMMON/SHARES/XUOO)  .  OEL(IOO)*  A  (lOOtlOO)  «NfMt  MN»NP1,NM1 

000003 

COMMON/EXPOPX  /  NEXOP1*  NeXOP2*  XeP1*  XeP2 

000003 

COMMON/  STiRXS/XsTR<  loo  )  *  XSSS<  l00)t  ODLL^OO) 

OOOOO3 

DO  10  j* 1 t  N 

0031 2q 

000005 

10 

xSSS(J)rX(J) 

003130 

mm 

DO  5o  J*l #N 

IF  tj.FG.1)  GO  To  20 

S8iUS 

000014 

I  T 

jM  1  *1-  1 

003160 

OOOO15 

X(jm*)«XSSS(jm1) 

003l70 

OOOOI7 

20 

x  ( j)  *XSS5  <  j)  ♦  Xppl 

0031 rO 

000022 

CALL  G^An*  S ( ]N) 

000023 

DO  3°  i-l.N 

003200 

0000?6 

30 

DDlL<I>*OEL(I>  „ 

00321 0 

000032 

X(J)«xSSS(j,.XEPl 

003220 

000035 

call  OHADlS(IN) 

000036 

DO  «o  I * J f N 

003240 

0000*1 

40 

*<0Di-1-<i>-del<i>>/<2**xepi> 

0032*>0 

000057 

50 

c0Nt,Nue 

003260 

000061 

XlN^XSSstN1 

003270 

000062 

RETURN 

0032B0 

000063 

ENO 

0032*0 

31 


subroutine  estims 


000002 

000002 

000002 

000002 

000002 


000002 

000002 

000003 

000005 

00000* 

000011 

000014 

0000*6 

000017 

000020 

000024 

000026 

000030 


000040 

000047 


8888H 

000066 

000070 

000071 


000075 

OOOIO* 

000106 

000110 

000113 

000116 

0001*5 

000127 

000127 

000133 


000140 

000144 


000155 

000156 

000157 


OCTOBER  1970 


C 
C 
C 

C  ESTIM  PERFORMS  The  COMPUTATIONS  TO  ESTIMATE  THE  LAGRANGE  MULTIPLIERS 
c  AND  MAKE  The  FIRST-  and  SECOND-ORDER  ESTIMATES  OF  The  final  solution 
c  of  The  problem. 

cOmmon/ShaRes/xiioO) ,  OEL(IOO),  A(ioo*ioo) ,n.m«  mn.npi.nmI 
COMMON  /EQALS/H.  Hi.  MZ 

common  /optnss/nti .nt2.nt3.nt4.nt5.nt6.nt7.nts.nt9.ntio 

common  /VALUe5/F,G,P0,r9IGMA,  hJ(200),  rhO 

COMMON/CRSTS/OELX ( lgol *  DeLXq ( 1 oo> '  RHOIN.RATIO*  EPSI*  THETAo* 

1  RSig1*  g1*  Xl(lOO),  x2(loO),  x3 (100)  •  xr2(100),  XR 1 ( 1 00 > .  pr 1 , 

2  PR2, pi.  F*.  Rj1  (200),  DOTT.  PGRAo ( 1 09 )  < 

3  PrEV3»ADEL*»  "TcTR,  NURINl,  "PhaSe 
C0MMON/CONPAX  /  NF 1  *  NF2*NF3 
CALL  STORES 


[ jr 


20 


.  J ♦  DI*G< 100) 
nSATIS 


Zio»RATIO**2 
Z9«RAT  10 

Zl«1.0/Z9*1.0/Z10 
Z2«Zl*l./Z9**3 
Z3»l ,/Z9**3 
Z4«Zl0*Z9 
j5»z9*«3 

46«1.0/( f  Zl 0-1.0) • (29-1.01 ) 
27«W29 

ze»i ./ (Z9-1  •) 

R0«1 «0/RHO 


IF  (NUMINI-2)  ,5n.B0M 
WRITE  (6,3J0) I 
PO«tPR2-Z**PRl*Z5*Pl)*. 

G«(RATl°*Gl-GRl)/(RATl8 

DO  20  I*1*N 

X  1 1 ) ■ (XRZ ( I ) -Z4*xr1 ( I ) *Z5*xl ( j ) ) *Z6 


10  CONTINUE 


5 


1.) 


NPsNPhASe 
NPHASE.4 
CALL  EVALUS 
nphase«np 
CALL  OUT PUX  (2) 

C  CHECK  TO  See  IF  estimates  have  CONVERGED 
GO  TO  <7n»30i7o>  »  NPHASE 
30  DO  50  j* 1 » M 

IF  iRJ(J))  *0,50,50 
40  IF  (THETaO‘rJ,JI ’  7o*50*5o 

5°  CONTINUE 

GO  TO  (7o»70*6Q> ♦  NT* 

60  C*LL  FINAlSINF3) 

70  CONTINUE  . . 

I  BO  WRITE  (6.340)  I  80  CONTINUE 

G«(RATIQ*GI-GRl)/ (RAtIO-1.) 

po«iz9*pi-p-i)*ze 


90 


DO  90  1*1. N 

X  ( I ) ■  (Z9*XI ( I ) -XR1 (I) )*ZB 

NP«NPhASe 

NphASe«* 


CALL  EVALUS 
NPHASE-NP 


003310 

003320 

003330 

003340 

88H18 


003470 

0034B0 

0034-0 

003500 

003510 

00352o 

003530 

003540 

003550 

003560 

003570 

00358o 


003600 

003610 


0036*o 

003650 


003670 

003690 

j;)?!; 

003720 

88!”« 

003750 

OO377O 

003B00 

003790 

003B10 

003820 

00383o 

003840 

003860 


32 


000161 

CALL  OUTpux  (?) 

OOSgflO 

C  CHECK  TO  SEE  IF  ESTIMATES  HAVE  CONVERGED 

00016? 

GO  TO  (140.100.1*0) .  NPHASE 

003890 

000171 

100 

00  1 2o  J.l.M 

00  3900 

000173 

IF  (RJ(J))  110,1?0,1?0 

003910 

000175 

110 

IF  (RJ(  j)  *thetao)  i*nM?o*i?o 

003920 

000200 

120 

CONTINUE 

OO393O 

000203 

130 

GO  TO  {  1 4 0  *  1 30  *  1 4 0 )  *  ^T4 

0039*0 

000212 

C*LL  FIN*LS(NF?) 

003960 

000214 

140 

CONTINUE 

000214 

1150 

WRITE  (6*350)  I  150  CONTINUE 

003?8o 

000220 

IF  <M)  1»0.1»0,1<>0 

000??? 

160 

D°  170  Jal,M 

003990 

888SI8 

128 

RJ( JJ.RHO^RJl ( J) 

IF-  (MZ)  21 0 f  c 1 0f  I9O 

88(8(8 

000?3? 

190 

DO  200  j*l*«Z 

00*0?0 

000234 

mKj.m*j 

00*030 

000236 

000^44 

2oo 

HJ (MNJ) *?.«Rji (MNJ) *RO 

00*0*0 

2io 

GO  TO  (220,f*0) .  NT2 

0°  230  I,1.N 

X  <I)*RH0'X1  <I> 

00*0^0 

0002*52 

220 

00*060 

000254 

230 

00*070 

888??? 

240 

CALL  OUTPUX  (2) 

CALL  REJECX 

000263 

250 

IF  (NUMINI-?)  ?80.300.?50 

00*100 

0002^6 

GO  TO  (?B0.310,?60) .  NT 7 

00*110 

C 

Sec°*0  orde«  M°VE  for  nexT  minimum 

00*120 

8881” 

l?8 

DO  270  I«1  •* 

DELX  (I)  “Z^X1  <i)-Z2*Xr!  ( I)  *Z3*xr?  ( 1 ) 

Pr2»pr1 

88^48 

000310 

200 

00*150 

00031? 

GR2,gr1 

OO4I6O 

OOO3I3 

PRl.pl 

004l?0 

000314 

GHI.gI 

OO4I0O 

000315 

DO  29°  I.l.N 

OO4I9O 

000^1 7 

XR2 (  1 ) *Xr1 (I) 

XR1 (I).Xl  (I) 

00*200 

0003?1 

290 

00*210 

0003?* 

RETURN 

00*220 

000325 

00033* 

300 

3 1  o 

GO  TO  (280.310.310).  NT7 

DO  3?0  I.l.N 

00*230 

00*2*0 

000336 

320 

DELX (I j .(XI (I)-XR1 (I ) )*Z7 

00*250 

0003*3 

C 

GO  TO  28() 

00*260 

00*270 

0003** 

330 

FORMAT  (/26H0  2ND  ORDER  ESTIMATES  ) 

00*280 

888?:: 

318 

FORMAT  </26Hn  1ST  ORDER  ESTIMATES  ) 

form4t  (/25H6  lagRange  multipliers  ) 

88*§o8 

000  3** 

ENO 

00*310 

000  002 
000002 
000002 
000002 
000002 


OOOOO4 

OOOOO5 

ooooo6 


000016 

000024 

000026 

888828 

OOOO4I 

0000^3 

0000*5 

000047 

000052 

000052 

88881? 

000057 

OOO06* 


000073 

000076 


OOOO75 

OOO16O 

ooolos 

000112 

000113 


SUBROUTINE  EVALUS 


C 

c  October  1970 

c 

C  IN  THE  NORMAL  phase  EVALU  CALLS  THE  USER’SUPPL IEO  ROUTINES  TO  EVALUATE 
c  the  objective  function  *no  the  constraint  functions  the  current 
c  POINT.  IN  THE  feasibility  phase  this  routine  PUTS  the  NEGATIVE  sum  of 
c  the  VIOcaTED  CONSTRAINTS  IN  LOCATION  F. 

cOMMON/SHARtS/X( 100) ,  OEL(100)f  a < 1 00  *  1 00> *n* m*  MN » NP 1 » NM 1 

£8mm8B  /^T^SS/NTyt1ftT2^NT3f  Nt4  '  NT5,  NT6f  NT7f  NTflt  NTy  f  Nti  o 
COMMON  /VALUES/F.G.POf RSIGMA,  RJ<200)»  RHO 

COMMON/CRSTS/DELXdOO) *  OELXOtlOO).  RHOIN. RATIO.  EPSI.  THETAO* 

1  §sIGl’  Gl*  XI ( 100)  .  X? ( 100) •  x3 ( 1 00 ) .  XR2 ( 100) .  XRl(lOO).  PR1. 

2  Pr2,P1,  FI,  RJ1 (200 ) ,  0°TT,  PgrAO(1oO),  OIAG(IOO), 

3  PREV3»AoELX*  ntctr*  NUMINI*  nphase*  nsatis 


c 

c 

c 

c 

c 

10 


H*0.0 

rsJGMA-°.° 

NSaTiS*2 

GO  TO  (10, 1°°, 190.200} .  NPHASE 

I'  E,6SSi®,LITT 

■3  GUESS 

•*  ALL  FUNCTIONS  ARe  TO  BE  EVALUATEO 
feasibility 

GO  TO  (20. *0) ,  NT2 


C  N0N-NEG4TiVieS  INClUOED 
20  0°  30  1*1, N 

IF  <X<I>>  260*260*30 
30  RSIGMA*RSIGMA-RHO*ALOG(X (I) ) 

IF  (M.eO,°)  go  TO  90 

oO  8°  J*1,M 

C*LL  rESTNX  (j,Rj(J)) 

IF  (RJM  j)  .LE.0.0)  GO  TO  50 
IF  (RJ(J) .GT.0,0)  G"  T°  60 
C  VIOLATION  OF  A  PREVIOUSLY  SaTISFIEO  CONSTRAINT 
GO  TO  260 

5o  IE  <RJ(J) .GT.o.O)  GO  TO  7o 

C  ALL  VIOLATEO  CONSTRAINTS  ADDEO  INTO  OBJECTIvE  FUNCTION 
F*F-RJ( J) 

GO  TO  8° 

60  RSlGMA*RSIGMA-RHO*ALOG(RJ(J) ) 

GO  TO  BO 

c  inoicates  satisfaction  of  constraint (iORmOre, 

70  nsat 1 s* 1 

RSigma»RSIGHA-RHO«ALOG<rJ(J)  ) 

80  CONTINUE 

90  C0nTiRue 

C  EQUALITIES  NOT  COMPUTEO  IN  FEAS.  PHASE 
PO.p.RSlGMA 
G*f-rhQ*fLOAT (M) 

IF  (NT2.EQ.1)  G»G-rHO«FLOAT(N) 

return 

c  regular  phase 

100  GO  TO  t 110’ 130) *  NT2 

C  NON  NEGATIVITIES  INCLUOEO 


00*330 

00*3*0 

00*350 

00*360 


00*370 


00*380 

00*3,0 


00*560 

00*510 

00*520 

00453o 

00*5*0 


00*570 

OOAKgO 

00459o 

00*600 

00*610 

miss 

004640 


00^66o 

004670 

004680 

004690 


00474° 

00*7*0 

004760 


00**90 


004800 

004810 

004820 

004830 

004840 

00485q 

004860 


34 


000121 

no 

DO  120  1 z 1 • N 

004890 

0P0123 

IF  (  X  ( I ) )  c60#?60#l2o 

00*^00 

000)25 

1?0 

RSIGmAsRsIGMA-ph0*AL0G( * 1 I  1  1 

00**10 

000135 

130 

IF  <M.E°#0)  G°  TO  15o 

004920 

mm 

DO  140  J*I*M 

CALL  RESTNX  (J.RJ(J)) 

004930 

00**50 

0001*2 

IF  ( R  J  (  J ) .LF.0.0)  GO  TO  ?60 

RSIGHA=RSigMa-»hO#aL0G(rJ{ J) ) 

000144 

004960 

000151 

HO 

CONTINUE 

004970 

C  EVALUATE  AND  ADO  IN  EQUALITY  CONSTRAINTS 

004980 

000154 

150 

CONTINUE 

004990 

000154 

CALL  RESTNX  (  o»F) 

0050I0 

OOOl^6 

IF  (m2)  l60'ieo,l6o 

000160 

160 

D°  170  1*1. MZ 

005020 

000162 

Js  1  ♦M 

005030 

000 1 £4 

C  ADO 

CALL  RESTNX  { j.Rj{ j) ) 

INTO  ThlRD  TERM  OF  P  FUNCTION 

ooSoSo 

000166 

HsH ♦ ( R J ( J) )**2 

005060 

000171 

170 

CONTINUE 

005070 

min 

oooi77 

I  00 

h*h/rho 

P°sRSIQMa*H 

P  0*p  4p  0 

8818^8 

OO5IOO 

000201 

G*2f *H-RhO*FLOaT (M) 

OO5IIO 

000205 

G  *  G  4  F 

005120 

000207 

if  (NT2tFQti)  g*g-rho*float (N) 

005130 

c  DUAL  VALUE 

005l*o 

000214 

Return 

005150 

C  GUESS  PHASE  NOT  COnEO 

005160 

000215 

190 

RETURN 

005IT0 

c-~ 

STRAIGHT  FUNCTION  EVALUATION  (  MATN^FEAS  ONLY) 

#81158 

000216 

200 

continue 

000216 
00021 7 

IF  (M.EQ.O)  GO  TO  220 

DO  21 0  1*1 ,M 

005200 

005210 

000221 

CALL  RESTNX  CXtRj(X) ) 

00523o 

000223 

210 

CONT 1NUE 

000226 

220 

CALL  REsT^X  (0,F) 

C  EQUALITY  CONSTRAINTS 

005250 

000230 

230 

IF  (NZ)  250.250.230 

005260 

000?32 

00  ?40  1*1, MZ 

00^270 

000234 

KZ«m*I 

0  05200 

888IJ3 

35?40 

CALL  RESTNX  (  KZ*KJ(KZ>> 
return 

005300 

c  CONSTRAINTS  violated  NOT  SO  BEFORE 

005310 

88*m 

260 

NS AT  I S*3 
p  0* 1 ® . £  35 

*©18 

0002*7 

return 

0053*0 

0002*7 

END 

005350 

35 


subroutine  feass 


000002 

000002 

000002 

ooooo? 


000002 

000003 

000011 

000012 

00001* 

000016 

8888*1 

00002* 

000032 

000033 

888835 

0000*1 

0000*3 

0000*5 

0000*7 

000051 

00005* 

000055 

0000*1 

000062 

00006* 

000066 

000067 

OOOO7O 

000072 

00007* 

000077 

000077 

000103 

00010* 

000105 

000105 

000105 

000105 


C 

C  OCTOBER  1970 

C 

C  FE»S  DETERMINES  WHETHER  THE  STARTING  POINT  IS  FEASABLE.  IF  IT  IS  NOT, 
C  FEAS  LOOKS  FOR  A  FEASABLE  ONE,  IF  NONE  EXISTS*  A  MESSAGE  IS  PRINTED 

c  and  control  returns  to  ma in, 

cOmmON/ShaRES/x (100) ,  OEL(IOO)*  A  ( X  00 • 100 ) ,N,M,  MNfNPlfNMl 
COMMON  /0PTNSS/NT1 • NT2 * NT3 *  NT* * NT5 • NT6 *  NT  7  * NT8 *NT 9 1 NT  1 0 
COMMON  /VALUES/F*G*Po*RSIGMA.  RJ<200>*  rho 

COMMON/CRSTS/OELXdOO)  t  OELXOdOO),  RH0IN*RATI0*  epsi*  ThETAO* 

1  rSIG1,  Gl,  XI (100) ,  X2 ( 100) ,  X 3 (100) •  XR2(100),  XR1(100),  PR1 

2  PR2*P1 *  FI*  RJ1 *200’ *  OOTT*  PGRAO<100>*  DIAG'IOO'* 

3  PREV3 * AOELX ,  NTcTR,  NUMINI*  NPHaSE*  nsatis 
NPH*SE*1 

G°  TO  ,lo,50) ,  NT 2 
10  NFIX'l 

00  30  1*1 * N 
TF  ( X  ( 1 ) )  20*2o*30 
20  *VlX«2 


X<I>»l.E-o5 
30  CONTINUE 

GO  TO  <5o**0>  *  NF I X 
*0  NPhASE“* 

CALL  EVALUS 

C  JUST  GET  ALL  CONSTRAINTS  EVALUATEO 
NPHASE*! 

WRITE  (6, I3U) J 

"c^lloUTpUxT^i 

IF  (M)  90,9**60 
DO  70  1*1  ,M 
IF  (RJ(I))  100,100,70 
70  CONTINUE 

80  CALL  TIMECS 

[WRITE  (6,1*qT1 


50 

*0 


6*0,0 

CALL  R£STnX  (  0»F ) 


CALL  OUT PUX  (2> 

90  RETURN 

100  CALL  BOOTS 

OO  110  I«1,M 
IF  (RJ( I) )  120,12o,110 
110  CONTINUE 

_ GO  TO  an 

I  I2n  WRITE  *6 ,1 50 n  120  CONTINUE 

c  to  indicate  to  main  to  st*rt  on  next  problem. 
nphase*5 

GO  TO  9° 


, 


c 

130  FORMAT  <1H0,2X,*8HMADE  VIOLATED  NON"NEGATI VI T IES  SLIGHTLY  POSITIVE 
1) 

1*0  FORMAT  <51h°*****THE  FEASIBLE  STARTING  POINT  TO  be  used  IS 
150  FORMAT  (3X,»9hTHIS  PROBLEM  POSSESSES  NO  FEASIBLE  STARTING  POINT,  W 
1  ILL  LOOK  FOR  DATA  FOR  NEXT  PROBLEM,  > 

EN0 


005370 

0053B0 

005390 

005*00 

00s*10 

00=*2o 


005500 

005510 

005520 

OO553O 

0055*0 

005550 

88H78 

005580 

00559() 

005610 

005620 


005650 

00566q 

005670 

005680 


00*710 


OO574O 

OOsjjO 

O0*77o 

00*78q 

005790 


005810 

OO502O 

005830 


005840 

005850 


005860 

OO587O 


005880 

00*890 


005900 


36 


SUBROUTINE  FlNALS<N?) 


000003 

000003 

000003 

000003 


000003 


0000?! 

000025 

000030 

OOOOaO 

000043 

0000*6 

000047 

000055 

000056 

8888£? 

000062 


OCTOBER  1970 

final  contains  the  tests  useo  to  oetfrmine  whether  A  point  satisfies 

The  FINAL  CONVERGENCE  CrItERION  CHOOSFN  TO  UETfRMINE  IF  The  NLP 

pkoblem  has  been  solved* 

n2  Set  FOUaL  TO  1  IF  CONVERGENCE  CRlTEplON  IS  SaTISFIEO* 
n2  SeT  eOUaL  TO  2  OTherwISe* 

COMMON/SHapES/X(100)  f  _OEL(100)  •  A|100.*  i°°i*N*M*  MN*NP1*Nm1 


10 

?0 

30 

40 


*  OFLXq (100*  *  RH0lN*RAT10*  tPSI*  THETAo* 
)•  X?(10°)t  x3(l°0)t  Xfi2<I°0).  xpl(100,t 


COMMON  /OPTNSS/NT1 .NT?*NT3.NT4«NT5.NTb»NT7tNT«tNTytNT10 
C°MmON  /VALuES/F,6,t>0.RSlGMA,  RJ(300),  RH° 
COMMON/CRSTS/OFLX(  ‘ 

1  KSlfi1!  G1*  Xl « 

2  po2,Pl.  F^t  R  jl < 2001  *  OOTTt 

3  rPEV3,AOELx,  rTcTR,  NUMlNj, 

GO  TO  llo^O'So1  *  NT5 
EPSIL*4BS(F/G-) .) 

IF  <FPS1L-T"ETA°)  50,50,70 
IF  (ABS(RSigmA)-TheTAO)  50t50.70 
IF  (NUMlNl_l)  50,40,40 
PEST=PR1- <PHI-P0,/ <1.-1 ,/SOHT  <RAT10>  > 


PR1 


PGHAO  <  1 0  0 )  «  OlAfidOOlt 

^PhaSF,  "satis 


EPSIL*ABS(PtST/G-l.) 

IF  <ePS1l-ThETaO)  50,70,70 

5o  n2=1 

G°  TO  (BotGO) ,  nT6 
(SO  CALL  PUNCHS 

GO  TO  80 
70  N2«? 

8°  RETURN 

END 


005930 

005930 


005940 

005950 

0059*0 

OO597O 

00591,0 

OO599O 


006070 

881858 

006100 

006110 

006130 

00*130 

006140 

00<>l5o 

006160 

006180 

88**88 

006310 
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000003 

000003 

000O03 

000003 

000003 


000003 
00001  I 


000024 

000026 


000031 

000037 

00004J 

000044 

000O51 

0000?6 

000061 

000O6I 

000063 

OOOO*4 

000065 

OOOO7I 

000072 

000073 

000076 

000101 

mm 

000 1 0^ 

88819f 

000113 

000120 

000122 

00012* 

000125 

000133 

000136 

000141 

000U4 

000146 

000155 


10 

20 

30 

40 


SUBROUTINE  GRABS (IS) 

OCTOBER  IB70 

GRAD  COMPUTES  THE  GRAOIENT  OF  THE  PENALTY  FUNCTION  ANU  THE  OUTER 
PRODUCT  FACTORS  OF  THE  MATRIX  OF  SECOND  PARTULS  OF  P. 

IF  <IS=1>  accum.  matrix  of  ?nr  paRtials  if<is*2>  0ont 

COMM0N/ShARES/X (10°)  *  DEL(100)«  A('00*100) *NtM*  MN*  NP  * t  NM  * 

common  /eqals/h*  Hit  mZ 

cQMMQN  /OPTNSS/NT ltNT2tNT3tNT4*NT5tNT6tNTjtNTB*NT9tNT10 

cOmmON  /VALUES/FtG,POtRSlGMAt  RJ<200),  rhO 

rOMMON'CRSTS^OELXIIQn) •  OFLAnnoo*’  RHOIN'RATIO*  EPSI*  ThfTAo* 

1  RSTg1’  g*’  x  1 C  *00) .  x2(l00),  x3(l00),  xr2(‘00),  xrI(IOO),  ppl, 

2  PFK.pl t  Fit  Rjl (2oo) t  ‘  ‘  ' 

3  prev3.a0ela,  nTCTR.  Nl 
GO  TO  (l0t30)  •  IS 
f>°  20  1,1, N 


OOTTt  pGR*D ( 1 oO ) t  0l*G(l00>* 
NUM IN  1 1  NPH asf •  NS AT I s 


DO  20  J«1*I 
A  < I • J>  =°  . 

00  0  I«1*N 

OELxO ( I ) ,0 ( 


c  this  section  works  correctly  in  feasibility  phaSf  as  well  as  normal  ph 

GO  TO  (SOtRU) t  NT2 
00  70  1*1 tN 
OELxO(!)a-RH°/x(I) 

GO  TO  <6o,7o) •  IS 
A (I . I)» ( -OEL X 0 ( I ) /X  (I) ) 


6° 

'0 

80 


IF  <M.LE#0)  GO  TO  180 

00  I70  k*i*m 
call  GRAoIS(K) 

IF  (RJ<K) .°.°)  GO  TO  110 

C  ALL  VIOLATEO  CONSTRAINT  GRAOS  AOOEO  TO  OBJ.  FUNCTION 
DO  100  I «  1 1 N 
IF  (OEL  < I >  >  90*  100.9o 
90  OELXO ( I ) «OEL XO  < I ) -OEL  < I ) 

100  CONTINUE 
GO  TO  170 
110  TT*RhO/Rj<K) 

00  l6°  I***N 


C  IF 
120 

130 


IF  (DEL<I>>  l20«lfe0#l2n 

oEL(it,o  skip  all  The  fallowing  c0mPutatiOn  inv°lvin< 

T*TT40EL(I) 

DELX°(ll«DELX°(I)-T 
GO  TO  {*30,460, f  JS 


G  «  BY  OEL ( I ) 


T"W° 


DO  1SU  Jj*!tl  4 
IF  ( OEL ( JJ) )  l40tl50.140 
140  A(t,JJ)«A(I,JJ) ♦T^OEL (JJ) 
150  CONTINUE 
160  CONTINUE 
170  CONTINUE 

C  equality  changes  FOR  GRAO 

180  IF  (MZ  «LE • 0 )  GO  TO  250 

GO  TO  (250,l9°t250) f  NPhASE 
I9o  RQ*2#/RHO 


006230 

006240 

006250 

006260 

006270 

0062B0 


00637q 

006300 


0064lO 

006420 

006430 


006440 

oo645o 

006460 

006470 

006400 

006490 

006500 

006510 

006520 

006540 

006550 

006560 

006570 

006500 

006590 

006600 

881118 

00663o 

006640 


006670 

00^^00 

006690 

006700 

006710 

006720 

006730 

006740 

006750 

006760 

00677® 

OO«>70O 
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000157 

000161 

mi:i 

000*67 

000171 

000172 

000175 

000?03 

000206 


000?l6 

000221 

000224 

000226 

000?J4 

000?36 

000246 


000255 

000?5? 


000263 


000264 


000275 


888??? 

000102 


000307 

000307 


HO  240  J=1 *^Z 
K  =  M*  j 

CALL  GRAplS(K) 

TTxRQ*R J ( K ) 

00  ?30  I r 1  * N 

IF  ( DEL ( I ) .EO.o.O)  GO  T n  ?30 
OELXO <1 ) =DELX0  <1 ) ♦  DEL  U 1 #TT 
GO  TO  (200,230) ,  jS 
?00  T*RO*DEL ( I ) 

D°  220  Jjrlil 
IF  (DEL  ( J  J )  }  210  *220  *210 
210  A(jtjj) *a ( 1 1 J J ) ♦T*PEL  <JJ) 

22o  continue 

230  C0nTiNue 

24o  CONTINUE 

250  GO  TO  (260*290)  *  IS 

?60  00  270  1  =  1. N 

27n  diag(i)*a(i*i) 

260  G°  TO  (290,330, 290 ) ,  NPHASE 
C  LEAVES  NEGATIVE  gradient  in  delp 
290  DO  300  1*1  * N 
300  OELXO(I)*-DtLXn(D 

310  ADELx*0. 

no  320  I * 1 • N 

320  AUEL X* ADEL X ♦OELXo  U) **2 

ADELX*SQRT (AOElX) 
return 

330  CALL  GRA01S<0) 

no  340  1 =  1 n 

340  OELXO  (I)  s-Dt-LXO  ( I )  -DEL  (J  ) 

C  LtAVES  The  nfg#  g«ao  °f  p  I*  nELxo 
go  TO  3 1 0 
END 


006790 

006800 

006820 
00683q 
006840 
006850 
OOftgfcO 
00^87  0 
0068P0 
QQ6890 

oo69oo 

006910 
006920 
006930 
006940 
006950 
OO6960 
006970 
0069R0 
006990 
0 0 '00 0 
007010 
007020 
007030 
007040 
OO7O5O 

00707  0 
007o9o 
007090 
007100 
007110 
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SUBROUTINE  INVERX(NSME) 

OCTOBER  1970 

INVERS  SOLVES  THE  SET  OF  EQUATIONS  FOR  THt  MOVE-VECTOR  USING  THE 
CR0UT  PROCEDURE-  If  The  Matrix  is  NOT  POSITIVE  DEFlNllE*  A  DIFFERENT 
METhOU  is  used. 


c 
c 
c 
c 
c 
c 
c 

PERFORMING  a  L"U  DECOMPOSITION  OF  THE  MATRIX  A,  TAKING  ADVANAGE  OF 

c**##the  symmetry  of  the  a  matrix. 

c#### if  a  non-positive  pivot  canoidaTe  is  generated  *  then  mccoRmicms 
c##*#PPOCEDURE  is  USeD(SeE  PP«  1 6 ?• i6g 


IN  FIACCO  Ano  MCCORMICK) - 


007130 

007140 

007150 

007160 

007170 

00T180 

007190 

007200 

007210 

0072*0 

00T23q 


C 

007240 

C****IF  NSME  »1  WORKING  WITH  A  NEW  A  MATRIX*  IF  NSME*  2  USING  PREVIOUS 

007250 

C*#** 

c 

a  Matrix*  rut  have  a  new  right-hand-side- 

007260 
00 >27o 

C****N1NV  is  The  nUMreR  Of  NOn.POsitIvE  pIV°T  CANDIDATES  GENERATED. 

0072B0 

000003 

COMMON^ShARES'X <100* *  DEL  *  1 00 * '  A <1 00  *  1 00 1  »N*M*  MN»NPl*NMl 

000003 

COMMON  /OPTNSS/NT1.NT2*NT3.NT4«NT5iNT6iNT7.NTb*NT9*NTio 

000003 

COMMON/CRSTS/DELXdOO) *  OELXOdOO),  RHOINiRATiOi  EPSI«  TheTAOi 

1  «SIG1,  Gl,  Xl<100),  X  2  (100) ,  X3(100),  X«2|100J,  X*l<100),  *»Rl , 

2  PR2  *P1 *  FI*  RJ1 *200  *  DOTT*  PGRAD'lOO'*  DIAG'lOo’* 

3  PREVs.ADELX*  nTCTR,  numini*  nphase*  nsatis 

0  0 
0  0 
0  0 
0  0 
0  0 
Ul  Ui 

COMMON/FXPOKX  /  NEXOPlt  N£XOP2*  XEP*t  XEP2 

DIMENSION  B(100) 

007360 

000003 

equivalence  <b*delx> 

007370 

000003 

2o 

GO  TO (20,  170) ,  NSME 

OO738O 

00*390 

0000 1 1 

NINVbO 

000012 

IF  (A(l,l))  40,30,50 

007400 

000014 

30 

NINV*1 

007410 

000015 

GO  TO  70 

007420 

88881* 

58 

NINV»i 

At1*1)*1. /At1*1) 

88?J28 

000021 

DO  60  I*2*n 

0t)7450 

000022 

60 

A(lfI)*A(l,l)#A(lfl, 

007460 

000030 

7  0 

DO  160  J*2  *  N 

007470 

888811 

OOOO35 

JM 1  *  J- 1 

t*6. 

qO  qO  1*1, jM 1 

887408 

007500 

000036 

IF  ( A ( 1  *  j) )  8o,9o*80 

0075i0 

000041 

80 

TeT^A ( J *  I j  #A ( I , J ) 

007520 

000051 

90 

CONTINUE 

007530 

000054 

A  t J* J) «A  < J, J) -T 

007540 

007S50 

000060 

IF  ( A ( j  *  j )  )  110*100*120 

000064 

100 

NiNv.ninv#1 

007560 

000066 

GO  TO  170 

007570 

lit 

NINV*NINV*1 

At J* J) *  1 . /A  <  J*  J) 

007580 

007590 

000074 

IF  tj.EQ.N)  GO  TO  170 

0076oo 

000076 

jpi«ja 

007610 

81818! 

00  150  L«JP1*N 

T*°. 

881118 

000102 

DO  UO  x*l ,  jmI 

007640 

000104 

IF  < A  1 1  *  j ) )  130*140*130 

007650 

000107 

130 

i»T4A(L,I)*A(I,J) 

007660 

000117 

140 

CONTINUE 

007670 

000122 

A(L*J)«A(L*J)-T 

007680 

40 


000127 

150 

007690 

000 1 36 

CONTINUE 

continue 

007700 

000140 

160 

007710 

00014? 

170 

CONTINUE 

007720 

00014? 

IK  <  N I NV )  1 B  0  *  1 00  *  2q0 

007730 

oooi44 

1«0 

B 1 1 » 

00  774  (j 

000146 

0°  ?10  J»?,N 

007750 

000147 

T*0  * 

007760 

OOOI5O 

JMlaj-l 

OO777O 

OOOI5? 

f)0  ?°0  I  X  1  «  JM  1 

007780 

0001^3 

IF  ( A  <  J • T ) )  l^Ot?00.l90 

007790 

000156 

190 

TbT^A ( Jf I )*B ( I) 

007000 

00016* 

200 

CONTINUE 

007010 

000167 

0001>4 

?10 

B<J>»<B<J>-T)*a<J.J> 

CONTINUE 

0°  ?40  I.1.NM1 

007820 

oo70"o 

000177 

007040 

000200 

NMKsN-I 

007050 

000202 

00  ?3°  J*1.I 

oo7B6o 

000?03 

L*NP1-J 

OO707O 

000?05 

IF  (A (NMK,L) )  ??0,?30,??0 

0078R0 

88811® 

158 

B (NMK) *9 (NMK) -A (NMK*L) #fl ( L ) 
CONTINUE 

QQ709q 

OO79OO 

000??1 

240 

CONTINUE 

007910 

000??3 

250 

GO  TO  <2B0*260>  *  NT3 

007920 

000?3l 

WHITE  <S»4^0> 

260  CONTINUE 

000?35 

?7o 

W«iTf  (6,*?o)  (DELXO ( I ) «l«l*N) 

000?50 

wHjTe  <6«440> 

270  CONTINUE 

000?54 

white  <6.4?0>  (OFLX(I) «I»ltN) 

000^67 

cBo 

RfTURn 

C^MPyTE  ORThOqON^l  MOV£ 

007970 

c— 

007980 

000270 

290 

continue 

007990 

0  0  0  27  Q 
000*7^ 

DO  350  1  1 s 1 • N 

000000 

I=N-H*1 

000010 

000276 

IF  ( A  ( I  ,  I )  )  31  0,300,320 

000020 

000102 

000304 

300 

H  ^ I } *0  •  0 

GO  TO  350 

000030 
0  Og04  0 

888581 

310 

B  ( I )  ■  1  #  Q 

GO  TO  330 

881818 

8881!! 

118 

B (I ) »0  •  0 

I  H 1  *  I  ♦  1 

IF  < IR1 *GTxN>  go  to  350 

DO  340  j»iP  tN 

881818 

000313 

000316 

340 

0  0  gO  qO 
OOglOO 

000317 

c0Nt  iNy^ 

000110 

000331 

350 

000120 

000334 

GO  TO  360 

000130 

c--  check  maybe  do  dief  for  p.s.d. 

88iu8 

OOO334 

360 

ZC?*°.° 

OOO335 

QO  370  1*1, N 

OOglgO 

000  337 

37o 

ZC?*ZC?*OELx° ( I) *B( I) 

006170 

000344 

IF  1 7C2 )  300,400,400 

000100 

000345 

300 

DO  390  I  * 1  * N 

000190 

000347 

390 

B  ( I ) s-B ( I ) 

000200 

000353 

400 

C  MCC 

WRITE  (6,450) 

zangwill  one  mod 

IF  <NEXCP2#NE,2)  GO  TO  350 

008??0 

000157 

000230 

000362 

DO  410  K*1 * N 

000240 

888181 

410 

B(K) *B  (K) ♦Dfc-LXo (K) 

GO  TO  ?50 

882118 

41 


000  37 1 
000371 
000371 


C 

420 

430 

440 

450 


FORMAT  (7E17.8) 

FFnDilT,<lH0,6X’ifHDEL  P  VfCTOR) 

rOo2^T  1H0,6x,24h§EC0nD  ORDER  MOVE  VECTOR) 

format  (iho,6x,i5hORth06onal  mo"e,  vtCT0R> 


00R270 

008280 

00fl2q0 

00®3q0 

008310 

008320 
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SUBROUTINE  OPTS 


000003 

000002 

000003 


00000? 

000003 

00000* 

88888? 

OOOoiO 


000103 

000105 

000111 

000113 

000113 


000011 

10 

ooool® 

000017 

20 

000021 

OOOO25 

30 

40 

8888?| 

c  MN 

OOOO3O 

OOOO3I 

000032 

000036 

50 

OO0O37 

OOOO46 

6® 

OOOO42 

000043 

70 

000050 

c  1  ' 

C  IF 

OOOOcl 

0000*0 

000062 

000063 

®0 

000064 

90 

8888 n 

00007Q 

OO0O73 

100 

no 

000075 

OOOlOl 

000102 

120 

MARCH  1971 

OPT  LOOKS  FOR  A  MINIMUM  ALONG  THE  SEARCH  VECTOR  USING  THE  GOLOEN 
SECTION  SEARCH  METHOD. 

COMMON/ShaRlS/XUOO)  ,  OEL(100*»  a<1°0»1°0> *N*M*  MN.NPltNMl 
COMMON  /VAi OeS/f»g*PO*RSigmA*  Rj(2001*  RmO 

c0mm0N/CrsTS/DELxi10OJ,  OELXO(IOO),  RHOlNtR*TlO,  EpSl,  thetao, 

1  RS1G1 *  61*  Xiiioo’*  X2  < 1 oo  *  *  xsfioo*’  XRStlOO*’  XR 1  *  1 00  *  *  PR1  • 
3  PR3.P1.  Fl*  R Jl (300 ) •  OOTT.  PGRAO<100).  01AG(100)» 

3  PReV3»Aqelx*  nTCTR.  NUmINI*  NPh*Se»  NSATjS 
kSU.I 
N405.1 
P3l»P0 

lSkl.i 

(jOTT»°. 

00  10  J*1»N 

dOTT»qOTT*oeLx ( J) *DELxO ( J) 

gQ  to  *0 

0°  30  I.l.N 
OELX  <1 1 b.DELX  < I > 

CONTINUE 
N40*“0 
MNbMN* I 

is  now  numb,  of  points  after  min  achIeveo 

NTcTRbNTtTR* 1 


•TRBNTrTR* 
50  J.I.N 


IS 


°9 

x2(i).x(I) 

PXi?Pn 
N  ^  0  1*0 

fg40 1  Bfg4®  1  ♦  1 
oO  7°  T  1 N 
X  ( i )  ■ ( i ) ♦OElX(i) 
CALL  EVALUS 


GO  TO  NSAT1S 

PX2«10,E35 
P0*l 0  tE35 
60  TO  10O 
CONTINUE 
PX2«Po 

IE  (PX1-PX2)  1 00# 1 00  f I5O 
IE  (N*0 1-2)  l3o*llOtA10 
DO  l2o  l«l#N 
Xl(I)-X(I) 

P1-PX2 
GO  TO  430 
ONLY  ONE  POINT 


so  far  computed 


150 


00  UO  I ■  1 1 N 

PpEV3.pXl 
GO  TO  lfio 
OO  160  1»1*N 


008340 

006350 

008360 

008370 

006380 


00®450 

006460 

008470 

oogcdo 

ooe^io 

00852o 

006530 

006540 

OOgSsO 


3857c 
°°85r° 

SS'lsS 

Q0861Q 

008620 

SS8632 

oo86*o 

00865° 

008860 

008680 

008690 

OO07OO 

00°' i0 

008720 

008730 

008740 

»o 

008770 

008760 

006790 

008800 

008810 

008820 

008830 

0088^0 

008850 

008860 

22111° 


000120 

00012* 

oooi?6 

000127 


000127 

OOOI3I 

000132 

000114 

000140 

000141 

000143 

000150 

000153 

00015* 


000166 
000 i 72 

00017* 

000203 

000205 

000213 

000?14 

000216 

OOO2I7 

000217 

000223 

milt 

0002P7 

000231 

888!?} 

000234 

000235 

000244 

000245 


000257 

000261 

OOO?63 

000264 

000267 

000271 

000276 

000301 

000301 

000306 


160  DtLX ( I) *1 .6lfl033Q9*OFLX (I) 

phfv3*px 1 
pxUpx2 
GO  TO  6° 

C  GOLDEN  SECTION  SEARCH  METHOD* 
r  R  vector  GOES  TO  X 1 ( I ) 

I7O  pOsl#^36 

N*  ® *sN*  O4 ♦ 1 
IB0  DO  I90  I=1*N 

l*o  xlmsxci) 

plA 

DO  200  I s 1 t N 

X ( I ) *. 38196601* (XI ( I ) -X3 ( I ) ) ♦ X3 ( I ) 

200  X2 ( I ) *  X ( 1 ) 

CALL  EVALUS 

GO  TO  <540*270*2101 •  NSATlS 

210  IF  (N404.LT. 30)  GO  TO  170 

211  continue 

c  THERE  is  NO  REFERENCE  TO  211*  THE  ABOVE  STATEMENT  IS  A  DUMMY  STATEMENT 
C  —  IT  iS  POSSIBLE  NO  FE*SlRLE  POINT  EXiST.  IF  NOT  TRY  MOVING  ON  DEL*0. 
C--  IF  IT  IS  NU j  POSSIBLE  TO  MOve  ON  DELX0  THE"  WE  MUST  BE  AT  A 
C--  SOLUTION  OF  NLP  PROBLEM. 

IF  (N404.GT. 100)  qO  TO  240 
220  D0  230  i31,n 

IF  <ABS(ABS(X3( j) /xl d) )-l.) .GT.l.E-7!  GO  TO  170 
230  C°NTlNUE 
240  GO  TO  <250  * 260 1 *  N4o5 

25q  N4  o5*2 

C--TRY  TO  MOVE  ON  GRaOIE"T. 

NTcTR*NTcTR-1 

MNxMN-1 

_ GO  TO  2° _ 

|26q  W^ITE  <6i~5BqT]  260  CONTINUE 

CALL  TimeCS 

call  OUTpUX ( 1 ) 

CALL  reject 
STOP  22042 
C 

27n  CONTINUE 

N4  04*0 

pxl*p° 

00  I*1 *N 

2B0  X (j) *0.3Bl^6601*(Xi (j).x2( j) ) *X2(i) 

CALL  EVAlUS 

GO  TO  <540*290.2201*  NSATlS 

290  llVxpJ 

300  N401aN40l»l 

IF  (N*01-25)  3*0 , 310,310 

310  kSW*2 

IF  (N401-*0)  3?0,*60,*60 
320  DO  330  1*1 *N 

IF  (ABS<X2(I)/X(I)-1.0> .GE.l.E-7)  00  TO  3*0 
330  CONTINUE 

GO  TO  460 

3*0  IF  (ABS(PX1/PX?-1 .) .UE.l.E-7)  GO  TO  *6o 
IF  (Pxl-PX^)  350 , 4  60  f  *00 

C  FHOM  LEFTORIGHT  x3(D (PREV3)x2(I) (PX1)X(I)PX2  X 1 ( I ) P 1 


00B900 

008910 

00fl920 

OO093O 

008940 

00fl9s0 

00fl960 

oofi97o 

00898° 

008990 

009000 

009010 

009020 

009030 

009050 

889878 

009071 

009080 

009090 

009100 

OO9IIO 

OO9I20 

0091 30 
0091*0 
009150 


009lp0 

0  O9 1 9o 
009200 


009250 

009260 

009270 


0093io 


009330 

88|?|8 

00,3fc0 

009370 

00938o 
009390 
009400 
0094 
009* 


1° 

?0 


009430 

009440 

009*50 

009*6o 


44 


00031 0 

350 

DO  360  IiI.N 

009470 

000312 

360 

XI <l>*x(l> 

0094R0 

000316 

c  THROW  AWAY  RIGHT  PaRT 

OO949O 

P1*PX2 

009500 

OOO3I7 

DO  37°  1=1. N 

0095 1 0 

C  POINTXPl  BECOME5*  Xp2 

00952o 

000321 

370 

X ( I)=. 381 96001*1 XI (I)-X3(I> >  *X3(I > 

009530 

c  temporarily  in  x  storage 

0095*0 

000330 

CALL  EVALUS 

000331 

GO  TO  <540«380* 170> •  NSATlS 

009560 

000340 

3fl0 

CONTINUE 

009570 

000340 

PXp*PXl 

c  switch  vectors  to  proper  position 

88115? 

000342 

Pxi *pn 

009600 

oo96io 

000343 

DO  390  I  *  I  *N 

000^44 

XX*X2(  I ) 

009620 

000346 

x2(I)mX(I) 

009630 

000347 

390 

x<i,sxx 

009640 

000353 

GO  TO  300 

009650 

C  LEFT  SIDE  TOSSED  AWAY 

009660 

c— 

CHANGES  for  nonunimodal  fn 

009670 

c— 

go  To  throw  a^a y  right  in  this  caSeinit  val  lt  fib  pt 

0096B0 

000353 

400 

IF  (PREV3-PX2)  350.350,*I0 

OO969O 

000^56 

*10 

DO  4?0  1=1. N 

009700 

000360 

x3(I|*x2(I) 

009710 

mm 

420 

X2  1 1 > =X  <1  > 

PRFV3»PXl 

88!!i8 

000  367 
000376 

430 

PXl=PX? 

00  **°  I.J.N 

009740 

OO975O 

000  372 

**0 

X(t)»0.38Iv6601*(xI ( I ) -X2 ( I ) > *X2(I) 

00976q 

000401 

CALL  EVALUS 

000402 

GO  TO  <540**50, 170> .  NSATlS 

009700 

00  04  1  1 

450 

CONTINUE 

009790 

ra:u 

PX?=Po 

GO  TO  300 

88||!8 

c  the  interior  points  now  give  equal  value  for  p,  compute  midpoint. 

009020 

000413 

460 

DO  47O  1  =  1, N 

OO903O 

0098*0 

000*15 

OELXO ( I ) «X  ( I ) 

000417 

X (I  Is (OELXO (I) *X2(I) ) *0.5 

009050 

000422 

470 

CONTINUE 

009060 

000424 

CALL  EVALUS 

000425 

GO  TO  (460. *90) .  KSW 

009000 

000433 

400 

IF  (ABS(P0/PX1-1.) ,gt.i.e-7>  go  to  5?o 

009090 

000441 

490 

GO  TO  (500,510) ,  ISw 

009900 

000447 

500 

IF  (Po.LT.P3I>  GO  TO  510 

009910 

OOO452 

ISW»2 

OO992O 

C  IF 

P-FUNcTlON  DlON.T  go  oOWN  TRY  NEG  VEcT. 

009930 

000*53 

GO  TQ  ?0 

ReTuhN 

0099*o 

000453 

510 

009950 

000454 

520 

DO  530  I*I*N 

009960 

88811* 

530 

X ( I ) =DELXo 1 1  > 

GO  TO  3=0 

888818 

c  are  we  now  in  feasibility  phase 

009990 

000ts$ 

000*6* 

5*0 

DO  550  I»I*M 

010000 

IF  ( R J ( I ) )  560|560f550 

oiooio 

000*66 

550 

CONTINUE 

n5aTiS,4 

0100?0 

000471 

010030 

000472 

RETURN 

0100*0 
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C - 

PROBLEM  has  hec.ome  FFASIRLE 

010050 

c 

P  -  FUNCTION  changes  if  a  constraint  become?  FEASIBLE 

0100^0 

0  0  O4  7  2 

56® 

MN  =  0 

OIOO7O 

OOO473 

5?0 

00  57°  I r  ^  f  M 

0100R0 

0004?5 

Rjl (I)sRj(I) 

0100^0 

000501 

C 

RETURN 

010100 

010110 

000501 

5A0 

FORMAT  (  floH  DPT  CANAT  FIND  A  FEASIrLF  POINT  t  THA  I  GIVES  A  LOWFR  V 
lALUE  OF  THE  P-FUNCHON.  ) 

010120 

010130 

OlOUO 

ooosoi 

END 
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SUBROUTINE  OUTPUX  (K) 

c 

010160 

c 

OCTOBER  1970 

010170 

c 

oioiflo 

C  OUTPUT  PRINTS  OUT  INFORMATION  ON  THE  RESULTS  OF  EACH  ITERATION  AND  THE 

C  SOLUTION  ESTIMATES  AND  THE  ESTIMATES  OF  THE  LAGRANGE  MULTIPLIERS 

010190 

010200 

000003 

000003 

C0mM0N/ShaRES/X(100> ,  DEL ( 100> *  A ( 100* 100> *N,M.  MN,  NP 1 » NM  , 

COMMON  /eOALS/h*  h'*  MZ 

cOmmON  /0PTNSS/NTl,NT2,NT3,NT*,NT5,NT6,NT7,NTR,NT9,NTl0 

000003 

000003 

COMMON  /VALUES/F,GtPOtRSlGMAf  RJ(200)t  RHO 

000003 

common/crsts/dflx < ioo> *  oelxo(ioo>»  rhoin, ratio,  epsi.  thetao* 

>  RSIG1*  G1'  xl (100) .  x?(l00>.  X  3 ( 1 0  0  >  *  XR2(100>*  XPl(lOO).  PRl. 

2  Pr2,p1,  FI,  rJ1(200),  D°TT,  rGrAO(100>,  DIAG(IOO), 

3  PREV3  ’ADELX  *  NTCTR’  NUMINI’  NPHASE*  NSATIS 

000003 

NZ»M*mZ 

010290 

000005 

GO  TO  (10,20, ,  K 

010300 

000013 

10 

WRITE  <6.60)  „  „ 

010310 

000017 

WRlTE  (6,70)  NTCTR,0°TT,RH0, ADELX, NpHASE 

010320 

| 888811 

20 

WRITE  (6  *8n)  F*P0*G»RSIGMa*H 

WRITE  (6,9°)  (X(I>*I»**N) 

818928 

000066 

WRiTf  (6,110) 

6°  TO  (3o*40) ,  NT2 

010350 

000072 

010360 

000101 

30 

WRITE  *6, 120' 

010370 

000105 

WRITE  (6.100)  (RJ ( 1 ) , I *1 »NZ  > 

01 0380 

0001?0 

60  TO  *0 

WrITe  (6,100,  (RJ(I, .I.l.^Z) 

0103r0 

000122 

40 

010400 

000135 

50 

RETURN 

010410 

C 

010420 

000136 

60 

010430 

000136 

70 

FORMAT  (iQX*6HP0INT»l4*6Xf 6H  D0TT*El5«7 »6X » 4HRH0«Ei5 t7 »6X * i (jHMAGNl 
1TUDE-E15. ' *?X.6HPhASE«12> 

010440 

010*50 

000136 

000l36 

eo 

9o 

FORMAT  (flXf2HF.El5t7,5Xf2Hp.El5.7f5Xf2HG.El5#7*5Xt7HPSlGMA«El5f7t5 
1X*2HH-E15.7  - 

010460 

miss: 

FORMAT  (6X#25hTHE  CURRENT  VALUE  OF  X  1S/(6e2°.7>> 

000>  36 
000136 

—•  O 
O  O 

FORMAT 

F^RMaT  i 6X * 21 HTHE  C0nSTRAInT  VALUES) 

0l04<)0 

010500 

mils 

120 

FORMAT  '28X  * 34HN0T  INCLUDING  THE  NON-NEGAT I V I T IES } 

END 

010510 

SUBROUTINE  OUTPUX  (K) 
9999  CONTINUE 
RETURN 
END 
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subroutine  pevalx 


000002 

000002 

000002 

00000? 

000002 


000002 
0000  03 

000004 

000012 

000014 

000024 

000033 

000034 

OOOO36 

000040 

000042 

OOOO47 

888813 

0  OOO56 
000056 
0000*0 
000067 
000071 

sswi 
000101 
000103 
000105 
OOO1 1 1 
000113 
000120 
000121 


C 

C  OCTOBER  1970 

C 

C  PEVALU  COMPUTES  THe  VALUE  OF  THE  PENALTY  FUNCTION  ANO  THE  VALUE  OF  THE 
C  DUAL  USING  PREVIOUSLY  COMPUTEO  VALUES  FOR  F *  *NQ  RJ. 

CQmMQN/Sh*RES/X< 100) *  OELdOO)»  * ( 100* 100) *N*M*  MNiNPliNMl 
c0mm0n  /eQaLS/h.  hi,  mz 

COMMON  ^ORYNSS^NTi *NT2»NT3*NTa»NT5*NT6«NT7,NT8,NT9*NTio 
COMMON  /VALUES/F,G,Po«PS1GMa.  R J ( 200 ) *  RHO 

COMMON/CrSTS/OELX ( 100) *  OELXO(lOO),  RH01NiR*TlO»  EPSl,  ThET*0i 

1  RSlGl,  Gl,  Xl(100),  X2 ( 1 00 )  ,  X  3  <100) ,  xRZllOOl,  xRl(lOO),  pRl, 

2  PR2*P1 *  FI*  RJ1 *200' ’  OOTT*  PGRAO<lo6’*  OIAG'IOO'* 

3  PREV3,AOELX,  NTCTR,  NUM1N1,  NPHASE.  NSAT1S 
H*0  •  0 

RSigmA«0,0 

c  nonnegs  if  Included  are  aodeo  to  p— are  pos,  in  all  phases 

GO  TO  (10,30),  NT2 
10  00  20  1*1 *N 

20  RSigMa,RsighA.RhO«al°G(X(1) ) 

30  GO  TO  <*o *50 *150) *  NPHASE 

c  objective  function  -  sigma  viol,  consts. 

*0  F*0 .0 

5°  if  (M)  loo, 100, 60 

60  00  9o  j«1,m 

IF  <RJ(J))  80,80,70 
70  RS1GMA»RS1GMA-RH0*AL0G*RJ( j> > 

GO  TO  9° 

80  F«F-RJ(J) 

9°  continue 

C  EQUALITIES  not  aooeo  in  FEAS.  phase 
100  CONTINUE 

IF  (MZ)  1*0*1*0*110 
llo  GO  TO  (1*0* 1^0* 1*0) ,  NPh*SE 

120  00  130  1«1,MZ 

K.M.l 

130  H«H*RJ(K)**2 

H-H/RHO 

1*°  hS«h*RS1gMA 

p°«f*hs 

HMS-2.«H-RhO«FLOaT(M) 

1FF(^T2,eQ,1)  G«G.RH°«FL°aT(N) 

1 50  RETURN 

end 


0105*0 

010550 

010560 

010570 

0105R0 


Ol067o 

010680 

010690 

OIO7OO 

OloTlO 

010720 

818^8 

OIO75O 

010760 

Olo77o 

010780 

010790 

010e00 

818^8 

010830 

0l0g40 

010850 

Olo86<) 

0l0870 

010880 

818*38 

OIO9IO 
0109?0 
010930 
0 109*0 
010950 
OlO’OO 
010970 
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SUBROUTINE  PUNCHS 


000002 

000002 

000002 

000002 

000002 


000002 

000002 

000004 

000027 

000030 

000043 

0000*3 

OOOlO3 

000113 

000114 

888M 

0001U 


c 

C  OCTOBER  1970 

C 

C  THIS  SUBROUTINE  PUNCHES  THE  STOPPING  POINTS  ANO  ASSOCIATED  PARAMETERS 
C  SO  THAT  ANOTHER  RUN  MAY  RE  MADE  STARTING  WHERE  THE  CURRENT  ONE 
C  STOPPEO 

C  THIS  ROUTINE  IS  CALLED  IF  NT6»o. 

cOmmon/Shares/x ( «  deli100)*  a(I°o*IO°) #N*M*  MN *NP 1 *NM 1 
COMMON  / EO ALS/h *  Hi*  M7 

COMMON  /OPTNSS/NTlfNT2*NT3,NT4*NT5.NT6*NT7fNT0fNT9fNTlO 
COMMON  /VALUES/F.G^O.RSIgMA,  Rj  (200,  *  RhO 

COMMON/CRSTS/DFLX  lioo' *  OELXO(100)*  RHOIN*RATIO*  EPSI*  THETAo* 

1  RSlGl*  Gl*  X  l  ( 100 )  •  X2  ( 1 00 )  *  X3  ( 1 00 )  *  XR2(100)*  XRK100)*  PRl* 

2  pr2.p1,  pit  Rjl (200) *  OOTT*  PGRAO(loO)*  OIAr(IOO)* 

3  PREV3.AQELXf  *TcTRf  NuMlNIf  nPHaSe,  nSATIS 
COMMDN'EXPOPX  /  NEXOPi*  NEXOP3*  XEPl*  XEP2 


T»a0,0 _  _  _  __ 

|  WRiTf  (7,1^)  rpSl  >RhO*  Theta0  *R~ATTO»TtMtNtMZ  | 

*mAx l  T S  SFT  Tn  * n •  SFTniunS 


C  TmmAX.jS  set  To  *0*  seconds 
NT  1  ,3 _ _ 

f^RlTE  (7,2qT  (X(I) tI»l.N) | 
C  SET  RHO  OPTION  SO  THIS  VALUE 


OF  RHO  WILL  BE  USE  FOR  THE  RESTART , 


WRITE  (7,3ft)  NTl tNT2tNT3*NT4tNT5tNT6tNT7tNT8tNT9*NTl0 
WRITE  C  7 , 2q )  XpPl*  XeP? 

WRiTf  (7,30)  NfXQP 1 »nfxqp2 

retu«n 


c 

10 


§8 


FORMAT  (5El2*5f 314) 
FORMAT  16E12.5) 
FORMAT  C  * 0 1 7 ) 

END 


010990 

011000 

011010 

011020 

011030 

011040 

011050 


011150 

0 1 1 1 7  0 
011190 

011200 


011240 

011250 

011260 

8H?r8 

0U2q0 
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SUBROUTINE  REJECX 


C 

011310 

C 

October  irto 

011320 

c 

011330 

c 

REJECT  RETURNS  The  STORED  VALUES  OF  THE  OBJECTIVE  FUNCTION*  THE 

011340 

c 

CONSTRAINT  FUNCTIONS  and  The  PENALTY  FUNCTION  TO  THEIN  NORMAL  LOCATION 

011350 

00000? 

COMMON/ShaReS/X(100) ,  DEL ( 100) t  A  ( 1 00 • 1 00 ) * N ,M*  MN*NPl*NMl 

00000? 

COMMON  /EOALS/H*  Hi*  MZ 

000002 

000003 

COMMON  /ValUES/F  *6*P0*RSI6mA*  Rj(200)*  RhO 
cOMM0N/CRSTS/OELxrl00) *  OELXO(IOO)*  RhOiH?RATiO,  EpSif  TheTAo* 

1  RSI  61  *  61*  Xi <ioo’ *  XZtlOO**  X3<l00)f  XR?  C 1 00 J  *  XR1(100,*  PRl* 

000002 

2  PR2.P1*  FI*  Rjl  (200 )  *  0OTT*  P6RaD<100)*  D1a6(100>* 

3  PREV3* AqELx*  NTCTR*  NUM1N1*  nPhASF*  NSATlS 

011*30 

DO  *0  1*1*N 

00000* 

10 

X  (1  )  *  X 1  (1  ) 

011440 

000010 

mmz*m*mz 

011450 

00001? 

qO  20  j* 1 f MMZ 

011460 

000013 

20 

RJ( J) *RJ1 C  j ) 
p0^1 

011*70 

000017 

011480 

0000?0 

RS16MA.RSIG1 

011490 

88SSI5 

6*Gi 

c*c  1 

8111?! 

0000?5 

~  '  « 

H*Hl 

01l5?0 

oooo?6 

return 

OII53O 

0000?7 

end 

0115*0 
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SUBROUTINE  RHOCOX 


000002 

000002 

000002 

000002 


000002 

88881? 

000016 

000020 

000021 

000022 


000024 

000025 


000034 

000036 

000040 

000042 


OOOO55 

OOOO56 

000057 

0000*2 

000067 

000071 

000072 

000073 

000074 

000075 


OOOO77 

000100 

000102 


000111 

000113 


000122 

00012* 

000130 

000136 


C 

C 

C 

C 

C 


£8 

30 

40 

50 

60 

C 


?0 


60 

9O 

100 

C  3 

no 

120 

c 

130 

c 

140 

150 

I60 

170 


OCTOBER  1970 


SUBROUTINE  TO  COMPUTE  INITIAL  RHO  VALUE 
CONTROLLEO  BY  COLt  7  ON  OPTION  CARO 

COMMON/ShaRES/X(100) #  DEL (100) «  a ( 100# 100) »N»Mt  MNtNPltNMl 
COMMON  /0PTNSS/NTl*NT2»NT3*NT4*NT5tNT6*NT7tNT6tNT9tNT10 
.Pn.RSiGMA.  rhO 


COMMON  /VALUeS/F  fG#PofRSlGMA#  WJ(200)# 

COMMON/CRSTS/OpLX ( 1 nn ) *  OELXn  ( 100 ^  #  RMOINtRATlO*  EP$1*  ThETAo* 

1  RSIG1*  G1#  Xl(*00),  x*(100)#  X3(100)#  xr2(I00)#  XrU'OO),  prI, 

2  PR2  *  pl  t  Fa#  Rj1(200)#  OOTT*  PGRAOUOO)*  01*0(100)# 

3  RrE V3 . AOELX »  «TcTR.  NuMiNj.  N^haSe.  nSaTiS 


j  BrEV3 • AOELX y  NTCTR#  NuMiNI# 
GO  TO  l110,5o,lOM9o)  *  NT  1 
RH0«RH0IN 
IF  (RHO)  30#3°f40 

8E?;a» 


NPARl*l 
RHO»j • 

2  means  rhd  which  minimizes  gradient  mag* 

CALL  GRaoS<0 
OO  7<)  i«*  «N 
pgR*d ( 1  *  »OELxO ( 1 ) 

RhO»2. 

CALL  GpApS (O 
OO  8°  I»**N 

OELXO ( 1 >  »OELXO ( I ) -PGR*0 ( 1 ) 
pGrAO(1).Pg“ao(1)-DELXO(1) 

GO  TO  (9o»l30>  *  NPAR1 
OOT1*0. 

OOT2-0. 

OO  100  I«1,N 

OOTl»OOTl*otLXO(l)«PgRAO(i) 

00T2»00T2«0ELX0(T) 

Rho»abs(oot1/oot2) 

G»  TO  20 

MEANS  COMPUTE  RHO  SO  AS  TO  MINIMIZE  OEL  P ''DDP'l . ' OEL  P 
NPaR2b1 
nPaR1»2 


USe  OF  and  Op  SUBROUTINE 
GO  TO  60 
RhO»i • 

assume  Sigma  term  is  cOnSio*  gRter  Than 
CALL  SECORX  (2) 

OO  1*0  l-l.N 
OELX(1,.pGrAO,1) 

CALL  INVERX  <11 
00  T5®  i.l,N 


Xl<D-OELX<l> 

OELX ( I ) sOELXO  ( 1 ) 

CALL  SECORX  (?) 

CAl 
00 

XR2(T)»0ELX<1) 

GO  TU  (170,200) ,  NPAR2 


ILL  SECORX  (2) 
ILL  INVERX  (1) 
3  U®  I»1,N 


OOT1*0* 


f  term 


0115*0 

011570 

0115BO 

011590 

011600 


0U6Bo 

nm 

OllTio 

011720 

011730 

81H18 

OII77O 

80 

oneio 

011620 

011630 

811*58 
011060 
01 1070 
0118*0 
011090 
0n900 
011910 

811|!8 

011940 

OII95O 

011960 

8HJJ8 

012000 

012010 

012030 

012040 

012050 


0120g0 

012090 

012100 

012110 
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mis 

O0T2*0. 

DO  ’8®  1*1 . N 

0001*1 

DOTl=OOTl«Pl>RAn(I  )*X1  <p 

01?140 

0001** 

180 

u°T?,o°T?.nELxn(i)*xp2(l) 

012150 

mill 

RHO  =  SORT(ARS(DOTl/OOT?)  ) 

G0  TO  20 

Mills 

0001S6 

1*0 

c 

NPAWps 2 

rho  minimizes  orqer  MOvr 

sist;s 

Oft  0 1 S7 

GO  TO  120 

012200 

c 

USES  INTERNAL  $UR.  TO  COM  /OOP/-1  Of  ANO  /OOP/-  OR 

nm 

0001^0 

?00 

OOT'-O.O 

oooi6i 

OOT?«0,0 

01223® 

000162 

00  210  1*1 ,N 

0122*0 

000163 

OOTl-Xl (jl**2*nOTl 

Ol221o 

000166 

?10 

0°T?axl  m*AR2(I).D°T? 

012280 

mr,i 

RHO*ARS<DOTl/OOT2* 

GO  TO  2° 

212270 

0l22n0 

000175 

ENO 

0122R0 
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SURROUTInE  SECOBX (IS) 


000003 

000003 

000003 

000003 

000003 


000003 

00000? 

ooooo6 

000017 

000OP4 

0000?6 

000027 

000033 

000017 

OOOO*5 

000047 

000060 


8888$g 

oooo*3 

000065 


00007* 

000075 

000100 


ooono 

000113 

000116 


000121 

000123 

000132 

00013* 

000136 

0001*0 

000U1 

0001** 

0001*5 

88811? 

000152 

000160 


C 

C 

c 

c 

c 

c 

c- 


OcTOBER  1970 

The  matrix  OF  SECOND  PArTIALS  of  the  PENALTY 


SECORD  EVALUATES 
FUNCTION. 

(1)  MEANS  DONT  COMPUTE  GRAD.  OUTER  PROOUCTflN  SECOHD) 

COmmON/ShARES/X  (  ^0)  •  DEL(l00>*  A  (  ^  ®0*  1 00)  #  n#  M*  MNtNP^NM1 


012310 

012320 

811128 

012350 

012360 

012370 


COMMON  /eQAlS/h*  h1*  mZ 

COMMON  /OPTNSS/NT l*NT2*NT3*NT4*NT5fNT6*NT7 
cOmmON  /VALUES/F.G^O^Slc^A,  RJ  (200 )  9  rhO 

*NT8*NTv*NT10 

COmMON/CRSTS/DELX(1qq)*  DFLX  Q ( 1 00  >  *  RHOIN* 
1  RSl  6  1  •  6^*  xl(100)#  x2(l00)#  X  3  ( 1  0  0  )  # 

RATIO*  EPS1*  THETAn* 
Xr2(100>,  XR^ ( 18°) .  prJ, 

2  PR2tpl#  Flt  R j 1 ( 200 ) •  DOTT*  pGRAD(lOO)#  D 

3  PrEV3,ADELX,  NTcTr,  NuMINI,  nPHASE*  nSAT1 

1 *G ( 1 00 » . 

no  io  1 B i * n 

01246Q 

10 

no  io 

0 1 247ft 

A  (p  j)*0. 

Ol2Afto 

GO  TO  (2 30*20) •  IS 

012490 

cgrag 

i.  TERM  NOT  PHEV.  COMPUTED 

012500 

20 

00  lo  1 a  p  N 

DO  30  j*l*i 

mi 

A ( I #  J) *0  #  0 

012530 

30 

CONTINUE 

0l25*n 

*0 

GO  TO  (40*60) *  NT2 

00  50  J  B 1  *  N 

012550 

01?560 

50 

A(I,I)«rH°/X(1,**2 

012570 

60 

continue 

0125*0 

]^i*-E.°l  GO  TO  130 

DO  120  IN*A*M 

811188 

IF  (R J ( In) )  120.12o*7o 

0l2*l 0 

70 

CALL  GRADlS<lN) 

TT  «  RHO^RJ < IN' ##2 

DO  110  IalfN 

811128 

IF  < DEL  < 1 ) >  8o.ll0.8o 

012650 

80 

tmtt#deu(i, 

012660 

DO  100  J* 1 • 1 

012670 

BO 

IF  (DEL ( J) >  ^0*100*^0 

A ( p j) «A ( j  « j) *T*OEL ( J) 

0126*0 

01?<>90 

100 

CONTINUE 

012700 

1 1  0 

CONTINUE 

012710 

120 

CONTINUE 

012720 

C  EQUALITY  CONSTRAINTS 

012730 

130 

IF  (MZ)  210.210.1*0 

0127*0 

140 

GO  TO  (2l0.150,230) .  NPHASE 

0l2750 

150 

««=2./BhO 

012760 

DO  200  JJ*1«MZ 

012770 

INsM.JJ 

0127RO 

CALL  GRADlS(IN) 

DO  iR0  I»1.N 

012800 

IF  <  DEL • I ) )  l60.1^0.1ft0 

012810 

160 

T«RO*DEL(I) 

012820 

1*0 

DO  lflo  J*1’I 

IF  ( DFL  (  J)  ) 

«o 

a1i»j)=a<i,j)*t*oel(j) 

cORTjNue 

012850 

180 

012860 

53 


000163 

000166 

000170 

000172 

000176 


miti 

00021* 

000?l5 


888IS8 

00022* 

000227 


000233 

000236 


88BIJ? 


000262 

000270 

000270 

000273 

00027* 

000276 

mm 

0003 1 2 
0003l5 
000322 
000324 
000*27 
000335 
000340 
000343 
000346 
000346 
000355 


88811* 

000362 

000364 

88831? 

000373 

000376 

000377 

000401 

000*02 

000*05 

000415 

000420 

000423 


190  CONTINUE 
200  CONTINUE 
210  DO  220  I  *  1 ? N 


0I*6(I)»A(I,I> 

220  A(I*T)»0, 

c  Ready  n0w  for  matrIx  of  2n0  p*rtial5  Of  Restrains 
230  GO  TO  (240*510’520> *  NT10 
240  IF  (N.LE.6l  GO  TO  340 
oo  330  in»**m 
.  LORN*2 

c  C0n5THa1nT  ASSUmED  n°nlineaR 
CALL  MATRXX  ( IN ’LORN ) 

IF  (LORN.LT.2)  GO  TO  330 
IF  (RJ(IN)  .GT.  0.0,  G0  To  2g0 
00  26"  l«z,N 

IMl»I-l 

00  260  j.lflNl 

IF  <  A ( j  t  j  >  >  250f26of?50 
250  A(I|j)tA|I|J)wA(JfI) 

A(J»l)«o. 

26°  CONTINUE 

00  270  X«1»N 
01Ag(I)*DIa^<I>“a<I*I> 


270 

280 

290 

300 


3l0 

320 

330 

340 

350 
C  — 


All*I)«0#0 
G6  10  330 
Ta-RHO/Rj( IN> 

OO  300  1*2, N 

i»ui-1 

0°  300  J.1.IM1 

IF  ( A ( J 9 1 ) )  290  *300 ’290 

A<I*J>«A(I*J)*T*A(J«I> 

A ( J* 1 ) *0, 

CONTINUE 

OO  320  I.l.N 

IF  (All. I))  31o  *320  *3 1 0 

01AG(l>»OlAG(I)*T*A(x,i) 

CONTINUE 

CONTINUE 

GO  TO  (520*350.520)*  NPHA5E 
IF  (mZ.Eq.O)  GOTO  420 
EOUALlTy  SeconO  paRTIAL5  HERE 
IE  <I}I10.GE*2)  GO  TO  420 
OO  410  1 1*  *  *  HZ 
IN«M*II 
L°BN.2 


CALL  MATRXX  llN’LORN) 

IF  (LORN.LT.2)  GO  TO  410 
T*2.*RJ( iN)/RHO 
OO  380  1«2.N 
IMl.I.l 

OO  370  J*  1  *  I H l 
IF  (A(J*1))  360*370*360 
36o  A(I*J)«A(I,J)*T*A(J.I) 

A ( J, I)«0,0 
370  CONTINUE 
380  CONTINUE 


012870 

012880 

012890 

012900 

012910 

012920 

81?’?8 

012,50 

012960 

012970 


012,,0 

013000 

013010 

81383$ 

013040 

013050 

8118^8 


013080 

0l3090 

013100 

013110 

013120 

013130 

0l3l40 

013150 


*lil?8 

0l3l80 

0l3l90 

013200 

°13210 

Oiazzo 

ol323o 

013240 

013250 

013260 

mi 


013290 


mm 
0 1 3320 
013330 


013350 

013360 

0»337o 

013380 

013390 

013400 

013*10 

013420 

013430 

013440 
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00  04?5 
000*27 
000432 


M43 
000446 

000451 

000452 

000*54 

000460 


DO  400  1 3 1 « N 

IF  (4  ( 1 1 1 ) )  39o,40o,3Qo 
390  0 1  AG ( I ) sOl  AG ( I ) «T«a ( I , I ) 

„„  All* 1*30,0 
400  continue 
410  CONTINUE 

c  get  matrix  of  2nd  partuls  of  objective  function 

420  ILL*? 

CALL  MATRXXI  °.LLL) 

IF  (LLL-LT.c)  GO  to  49o 
D°  440  I s2  « N 


013450 
0  I 346o 
0  1 347 0 
0134BO 
0l34qd 

013500 
01 35 10 
013520 

0  I 354o 
013550 


mai 

IMI *  I- 1 

DO  440  Ja 1 ♦ IM 1 

811158 

ooo4*4 

430 

IF  <A<j,j))  430l440f430 

0l35Rp 

000467 

A( I t j) a* ( I t J) *A { Jf I) 

0 I 

000476 

440 

A  (  J  ,  I  )  3  A  { I  ,  J) 

013600 

nm 

4S0 

00  470  1*1 *N 

IF  ( A ( I « I ) )  450,460,450 

811118 

000*516 

Ail*I)*DIAG<I)*A(IM) 

G°  470 

0 1 3  63  q 

000532 

0 1 364n 

000523 

460 

A  < I »  I  > *0 I  AG ( I ) 

013650 

0005?7 

470 

CONTINUE 

013660 

000532 

480 

RETURN 

013670 

000533 

000*35 

490 

DO  500  1*1  * N 

A<pi)aDiAG(l) 

0I36B0 

0l36qn 

000541 

Du  500  J 3 1 ♦ N 

013700 

000542 

500 

A(PJ)*A(J<  I) 

0 1 37 1 0 

000556 

GO  TO  4fl 0 

013720 

000556 

0005fi5 

^o° 

G0  IQ  <520*350*350) *  NPHASE 

DO  53 0  1  a ^  ,  N 

im 

000567 

IMl.I-l 

013750 

888;;i 

530 

DO  530  J  a 1  •  I M 1 

A(JiT)aA(T, J) 

DO  540  I  ■ 1  *  N 

811558 

000*06 

0 I 37O0 

000607 

540 

A ( I , I ) 3D  I  AG ( I ) 

013790 

0006)6 

GO  TO  400 

013600 

0006)7 

END 

01 36 10 
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SUBROUTINE  SETS (TMMAX) 


888885 


ooono3 

888885 

oooni2 

000013 

000014 


c 

C  FEBUARY  1971 

c 

C  SET  STORES  ThE  TIME  AT  WHICH  THE  PROBLEM  IS  BEGUN 
COMMON  /TSWS/NS^W 
COMMON  /TMXS/TM0fExT*EXT9n 
C 

c  sec0nd  gives  jOr  Cpu  execution  time  i*  i/iooo  °e  a  sec°nd 
c 

CALL  SECONO  <TMO) 

EXTsTMMAX*TMO 

eXT9°*  TmO  ♦  °.90#TMMAX 

NSWWal 

Return 

ENO 


013030 

013040 

013050 

013860 


O1309o 

013900 

013910 

0139P0 

813128 

0 1 395(1 
013960 

013970 
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SUBROUTINE  STORES 


000002 

000002 

000002 

000002 


000010 

000012 

000013 

000017 

000020 

000022 

000023 

8888?! 

0000?7 


C 

C  OCTOBER  1970 

C 

C  STORE  STORES  THE  VALUES  OF  THE  CURRENT  POINT  AND  THE  ASSOCIATED 


c  Values  of  the  functions  in  a  temporary  area. 

COMMON/ Sh ARES/X ( 100 1  *  DEL<100)  *  A<l00,l00>*MfMt  MN  *  NP1 «  NMl 
CQMMQN  /EOALS/Ht  H 1  *  MZ 

c0mm0n  /VALUES/F, G,p0tRSlGMAt  RJ(200)t  RH° 

COMMON'CRSTS'DELX ( Iqq) *  DELXq (100*  *  RhOIN*RATIO»  LPSI*  THETAq 


HS] 


51  *  g1*  Xl(l00)f  x2(10d)#  x3(l°0)»  Xr2(100)*  XrM100). 


s IG  *  G1* 

2  pR2fpi,  Fl*  R  jl  ( 200 )  t  DOTTt  PfiRAnUOO).  OIAG< 

3  PREV3#ADELXf  NTcTR,  NumiNi#  nPHA$F.  *SaTI$ 


100)  * 


10 

?0 


DO  10  I  *  1 • N 
XMIJ-Xt!) 

dO  20  j* 1 f mMZ 

P1*P0 

Fl.F 

G1*G 

RS I Gl “RSI GMA 

return 

end 


PR1. 


013990 

0U000 


014030 


812B8 

01*130 

ouno 

01*150 

01*160 

01*170 

oi*ieo 

01*190 

81*198 

01*220 
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SUBROUTINE  TCHFCX 


000002 

ooooo? 

000002 

000002 

ooooo* 


000007 

000016 

000012 

000020 

000020 

00002* 

000026 

000027 


FEBUARY  1971 

TCHECK  CHECKS  THE  NUMBER  OF  SECONDS  THAT  HAVE  ELAPSED  SINCE  THE  START 
OF  THE  PROBLEM.  IF  THE  SOLUTION  IS  TAKING  LONGER  THAN  90  PER-CENT 
OF  The  ESTIMjTeD  MaXjMUM  time?  A  S«ITrH  IS  Set  TO  GIVE  MORE  OUTPUT. 
COMMON  /OPTNSS/NTl*NT2*NT‘**NT**NT5tNT6tNT7,NT8.NT'**NTlO 
COMMON  /TSWS/nSVW 
c0mmon  yTMxS/TMOtEXTfEXT9o 
CALL  SECOND  <SECS> 

IF  (SECS.LT  »EXT9n>  RETURN 

GETTING  cLOSe  TO  EXCEEDING  The  Time  limit  Set  OUTPUT  option  TO  give 
more  OUTPUT. 

NT3«2 


X«SFrS  -  TmO 

|writf(6»»»)~x1 

10  FORMAT  jme“*F9*3»8H  sEC0N0$  > 

IF(SecS  #GT#  ExT)  nSww«2 
CALL  outpux  <1> 

RETURN 

END 


0U240 

014250 

0142^0 

014270 

0142R0 

014290 


014330 

811318 

01*360 

01*370 

01*3*0 

01**00 

01**10 

01**30 

01***0 
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SUBROUTINE  TIMECS 


000002 

000002 


88883! 

OOOOOfc 

00001* 

000017 

000020 

000021 

000021 


C 

C  FEBUARV  1971 

C 

C  TIMEC  CHECKS  THE  NUMBER  OF  SECONDS  THAT  HAVE  ELAPSED  SINCE  THE  START 
C  OF  THE  PROBLEM.  IT  PRINTS  THIS  NUMBER.  IF  THE  SOLUTION  IS  TAKING 

c  timec  checks  the  numrEr  of  seconds  th4t  h*ve  e<-apsed  since  the  st*rt 

c  OF  The  PROBLEM,  IT  PRINTS  this  NUMBER.  IF  THE  SOLUTION  IS  TAKING 

c  longer  than  the  estimated  maximum  time,  a  switch  is  set  to  terminate 
c  the  «un. 

C OmmON  /TSwS/NSww 
COMMON  /TMXS/TMO.EXT.EXT90 
C 

C  SECOND  GIVES  JOB  CPU  EXECUTION  TIME  IN  1/1000  OF  A  SECOND 
C 

CALL  SECOND  (SECS) 


10 

C 

20 


X«SfcS-ThO 

WRtTf  (6.c°)  XI 

IF  (SECS. LT. EXT) 

NS*w,2 

ret1,rn 


GO  TO  10 


FORMAT  (6X,5HTIME»F9»3.8H  SECONDS) 


END 


014460 

014470 

014480 

014490 

014500 

OU49O 

014500 

OU5IO 

01*S2o 


014550 

014560 

OU57O 


0l*6l0 

014620 

014630 

014640 

014650 

014660 
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SUBROUTINE  XmqvF S 


000002 

00000? 


000002 

000002 


000002 

000012 

0000U 

000016 

0000?0 
000  0?4 


000025 

000o?6 

000027 

OOOO3O 

000032 

000033 

000034 


0  000  35 
00003^ 
000037 

88831? 

0000*0 

000062 

88885? 
000  072 
000073 
000076 


000102 


MARCH  1Q7] 


1, 


XMOVE  DETERMINES  the  VECTOR  along  WHICH  THE  SEARCH  for  a  minimum  is 
USING  opt. 

cDMMON/ShaRES/X<100) «  OFL (100) «  A  (  100*  100) , N,M*  MNtNPl.NMl 
COmmON/CRSTS/OFLX(IOO) *  OFLXO(iOO)*  RHOI Nt R AT 1 0*  tPSl*  THETA0» 

1  WSIGI,  Git  XI(IOO),  X2  ( 1 00 ) .  X3  ( 1 00 ) .  XR2ilO(Mf  XRI<100).  PR 

2  HR2  *P  I  *  Fl*  RJl<200)f  OOTT*  PGRAO^IoO^  OIAG<l00>f 

3  RRFV3 • aoelx  *  NTcTR*  NUM INI*  NPHaSF*  nsatis 

CQMMQN/F.XPOHX  /  NEXOPl.  NFXDP2.  XFP1.  XEp2 
cOmmON  y  ^VF^/Sj  G ( I  00 ) , YY (I  00) t XXX (100) .DELL (100) 

C  — NEXOP?  DETERMINES  how  move  is  to  be  made 
c  nexor2  *  1  usE  mooified  nekton  raPhSOn  methOo* 

c  *  2  MOolFlEO  NpWTON  RaPhSON  METhOd*  pUT  ADO  DE*-*0  TO 

r  orthogonal  move  vector  if  hessja^  is  inofmnite* 

C  ,3  USE  STEepeST  DESCFnt  METh°0. 

C  3  4  USE  MCCORMICKtS  MODIFICATION  OF  THE  FLF TCHER-PDWELL 

c  METHOD. 

GO  TO  (I0«l0*lfl0.30) «  NEXOP2 
c--MEWTON  -RAPh  WiTh  WHATEVER  MeThOO  is  in  INVERSE 
10  CALL  GrAdS(1) 

C--ONE  (I>  MEANS  ACCUMULATE  matrix  of  second  partial  OERIVATIVES 
CALL  SECORX  (I) 
ftD  20  !*1#N 

20  OELX(D*OELXO<1) 

CALL  INVFRX(l) 

c  IF  A  NONPOSITVF.  PIVOT  IS  ENcOUNTFREO  IN  INVERSE  AN  ATTEMPT  IS  MADE  TO 

c  compute  fl  vector  having  a  positive  dot  product  with  a  negative 
c  EIGENVECTOR  and  THf  NEGATIVE  of  del  P. 

CALL  stdres 
CALL  OPTS 
RETURN 


C- -F-P-D-MCC  move 
3°  CALL  GRAoS(2) 

c--mn  is  no*  Of  moves  for  this  value  of  pho 
IF  (NN.Ne.O)  6o  to  ?0 
40  IREP*0 

IT,0 

C— SET  INITIAL  GUESS  INVERS  MATRIX  OF  SECOND  PARTIAL  DERIVATIVES 

c—  use  partial  inverse  if  known 
do  50  1 « 1  •  N 
DO  ’0  j* I  *  N 
50  A ( I t  J ) aO  #  0 

DO  60  1*1 *N 
6®  AdU)*1,0 

TO  DO  *0  1  * 1  * N 

*0  DELX ( 1 ) *DEL XO  ( 1 ) 

IF  (IREP.GT.N)  GO  TO  40 
IF  <IT.EQ.°>  GO  TO  130 
00  90  l*I,N 
SlG<l)*X(l)-XXXCl) 

90  YY(1)*DELL(I)-DELX0(1) 

c— NEGATIVE  gradient  stored  and  computeo 

C--COMPUTF  HY 

0D  100  |B1vn 


014680 

014690 

0147oo 

014710 

014720 


81i!8I 

oi*en? 
01**0'’ 
01*80* 
01*805 
0I*80A 
0I*pI0 
0  1*8?0 

01*8*0 

oi*88° 

Ol*87o 

01*890 
01*900 
01*91 0 


01*9*0 

01*9S0 

OI497O 

01498° 

014990 

015000 

015010 

OI5O2O 

OI5O3O 

0l’040 

015050 

811*;* 

0I5080 

015090 

sum 

0l5l?0 

0l5l30 

0151*0 

811118 

015I70 


60 


000104 

000105 

000107 

000122 

000123 

0001?* 

000132 


000136 

0001*0 

000143 

888115 

000165 

000167 

OOOI7I 

OOOI74 

000176 

000177 

000201 

000214 

000215 

000^16 

000223 

000225 

88815! 

000233 

888531 

000236 

888^8 

000?*2 

000246 

000247 

000250 

000251 


DELX (11*0.0 
DO  100  J*!.N 


100  DELxd>»DELxd>*A(I*J>*YY(j) 
C— C°MpUTE  Y(S1G  -HY)  -1 
ZCON*o,0 
do  no  i*i,h 


110 


c— 


120 

C  — 
130 


ZCON*ZCOn*YY(I)«(SIG<1>-DeLX ( 1 > ) 

IF  (ZC0n,E°.0.0)  G°  T<5  130 

1REP*1REP*1 

ZC-'./ZCON 

update  h  matrix  using  mcc  formula  when  scalar  ne 

00  '2°  1* 1 tN 

T1*ZC* (SiG  < 1 >  *DELX  <  I  >  > 

0°  120  J,1,N 

A(1*J)*A(1*J) *T1* (*OELX ( J) *S1G ( J) ) 

A ( J» 1 ) *A ( I » J) 

STORE  CURRENT  POINT  ANd  CURRENT  GRAOIEnT  <NEG> 

0°  140  1,1,* 

xxx<i>*xd> 


UO  DELL  ( I )  *oELx0  ( 1 ) 

00  ISO  1*1, N 
OELx ( 1 >  *0  •  0 
0°  150  J,l,* 

150  DELX (I>*OELx  <1> ♦ A ( 1  * J> *OELxO «J» 
ZC1*0,0 
OO  ^6®  1*1 .N 

16o  ZCl»0ELxd)**2*Zcl 
ZC1.S0RT(ZC1) 

OO  170  1*1, N 

1 7o  DEUdl)*OELxd)/ZCl 
CALL  STORES 
CALL  OPTS 


1  T  » 1 T  *  1 
RETURN 

180  CONTINUE 
C  STEEPEST  DESCENT 
CALL  GRAdSIO 
DO  *9°  1  *  * , N 
l’O  OElx(I)*OEi XOll) 
CALL  STORES 
CALL  OPTS 
RETURN 
ENO 


015180 
OI5I9O 
01 52 00 
015210 
015220 
01g230 
01*2*0 
015250 

811I?8 

0152R0 

015290 

Ol53oo 

015310 


Ol53*o 

015350 

015360 


811328 

01539o 

015400 

015410 

015420 

015430 

0 15**0 
015450 
015460 
01547Q 


811188 

015520 

015530 

015550 

0l5560 


015590 

015600 
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*> 


User-Supplied  Subroutines  for  Inside  Programs 


63 


000002 

000002 

000003 


SUBROUTINE  REAOIx 
9999  CONTINUE 
RETURN 
END 


64 


\  fVP.  '  f  1  ■n  •-  1  -w  1  7T  r  F  w  '  ✓  IT? Ir 


000005 

000005 

000005 

000005 

000007 

000011 

000011 

000013 

000015 

000023 

000024 

000025 


000033 

000035 

000035 

000037 

mm 

8888H 

000057 

000063 

000064 

000064 

000065 


SUBROUTINE  RESTNX ( INfVAL) 

COMMON/PROB/ISP 

C0MM0N/IN/W(l00) 

C0MM0N/SHARtS/X(100) »0EL<100) fA(l00t10°) .N, M, MN,NP 1 iNM1 
^(TSP-1)  1o5°,1000t2000 
1000  f"," 

VAl*0» 

IF (IN) 1100tll00»1200 
1100  DO  n50  j»l»N 
11=0  VAl»vAl*w(j)*(X(j)'2,)**2 
G°  TO  9999 
1200  VAL*FN 

_  DO  1250  J*1»N 
12=0  VAl»vAl-x(j) 

GO  To  9999 
2000  FN»N 

val*o. 

IF <1^)2100,2100,2200 
2100  DO  2150  J*1 *N 

2l5o  VAL  ■  VAL  ♦  W(J)*»,5  *  <X(J)-2»)**2 
GO  TO  9999 
2200  VAL«FN 

D0  2250  jcl «N 
225(j  vAj_»vAl-x  ( j)  **2 
GO  TO  9999 
9999  CONTINUE 
RETURN 
END 


65 


000003 

000003 

000003 

000003 
000005 
00001 0 

000013 

000015 

000017 

0000?5 


0000?6 

000030 

000034 

000034 


000052 

000053 

000055 

000061 

000062 

000062 

000063 


SURROUTlNE  GRAfllS (IN) 

COMMON/PROR/ISP 

COMMON/lN/WdOO) 

COMHON/ShAR^S/X ('00) ,OEL (100> *A (l00*100) ♦N*MtMN»NR1 tNM1 
DO  500  J= 1  * N 

500  del(j)=9. 

IF(lSp-l) 1000,1000,2000 
1000  IF  (IN)  1100.1100»l200 
1100  D°  1150  Jsl»N 
1150  DEL ( J)  *  2.  W(J)#(X(J)-?.) 

GO  TO  9999 
1200  DO  1250  J=1»N 
1250  DEL(J,*"1» 

G°  TO  9999 

2000  IF(IN)2100*2100»2200 
2100  DO  2150  J*1*N 

215°  DEL(J)  *  2.  *  W ( J) *#.5  •  (X(j)-2.) 

GO  TO  9999 
2200  D°  2250  J  =  1»N 
2250  DEL ( J*  x-2.*X(J> 

GO  TO  9999 
9999  CONTINUE 
RETURN 

end 
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000005 

888881 

000005 

000007 

looo 

oooon 

1100 

000013 

000023 

1150 

000023 

000023 

1200 

000024 

000020 

2000 

2100 

000030 

000044 

2150 

000045 

2200 

000047 

22->0 

OOOO56 

000056 

000056 

000057 

9999 

SUBROUTINE  MATRXX  (IN*L) 


common/prob/isp 
COMMON/lN/W( 100) 

COMMON/SHARtS/X(l00> ♦OEL(100) *A (100,100) , N*M*MN» NP1  .NM1 
IF(ISP-I) l0^0»1000t2O00 
IE  ( IN)  H0o,noo»l200 
D°  1150  J*1.N 
A ( J  * J )  a  2,  W (j) 

GO  TO  9999 
CONTINUE 


GO  TO  9999 

IE (IN )2 100*2100* 2200 
DO  2l50  j»**N 
A(J,J)»  2,»W(J)*#.5 
GO  TO  9999 


DO  2250  J“1 

A ( j» j) »-2, 

GO  TO  "99 
CONTjNue 


*N 


RETURN 


end 
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USER-SUPPLIED  INFORMATION  CARDS  AND  RESULTS  PRINTING  OUT  ALL  POINTS 
OF  INSIDE  PROGRAMS 

It  is  useful  in  debugging  the  overall  program  and  in  increasing 
solution  efficiency  to  print  out  the  solutions  of  the  inside  programs, 
preferably  all  points.  The  starting  points  of  the  inside  programs 
should  be  feasible,  which  is  dependent  upon  the  starting  point  of 
the  outside  program.  The  starting  point  of  the  outside  program 
should  also  be  feasible.  The  solutions  of  the  inside  program  for 
the  starting  point  of  the  outside  program  should  be  examined.  Also, 
several  points  of  the  outside  program,  and  the  corresponding  inside 
program  solutions,  should  be  examined. 

The  SUB  computer  program,  leaving  in  the  statements  for  printing 
of  points  of  the  inside  program,  provides  the  opportunity  to  accom¬ 
plish  the  above,  when  used  with  the  user-supplied  information  cards 
given  here. 

User- supplied  information  cards  are:  parameter  card,  initial 
starting  point,  option  card,  tolerance  card  and  second  option  card. 
Three  sets  of  user-supplied  information  cards  are  necessary  for  the 
example  problem. 

The  outside  program  has,  on  the  parameter  card,  information  that 
EPSI  =  l.E-05,  RHOIN  =10.,  THETAO  =  l.E-05,  RATIO  =10.,  TMMAX  = 

360.,  M=2,  N=4,  MZ=0.  The  initial  starting  point  is  x^  = 

=  x^  =  x^  =  4.  The  option  card  indicates  that  all  points  of  the 
outside  program  are  to  be  printed  out  (Option  3=2),  and  that  there 
is  at  least  one  nonlinear  constraint  (Option  10  =  1).  The  tolerance 
card  indicates  that  for  numerical  differentiation  and  control  of 
P- function  minimization  the  tolerances  are  .001.  Finally,  the  second 
option  card  indicates  that  the  problem  is  to  be  solved  without  check¬ 
ing  first  derivatives  (Option  1=1),  and  that  steepest  descent  is 
used  to  minimize  the  P- function  (Option  2=3). 
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The  inside  program  information  cards  are  similar,  with  the 
following  exceptions.  Each  problem  has  one  constraint  (M  =  1). 

The  starting  points  are  vg  =  v2  =  v3  =  v4  =  • 5 •  The  first  option 
cards  indicate  that  all  subproblem  solutions  of  the  inside  programs 
are  to  be  printed  out  (Option  3=1).  The  first  option  cards  also 
indicate  that  the  first  inside  program  has  all  linear  constraints 
(Option  10  =  2)  and  the  second  inside  program  has  at  least  one  non¬ 
linear  constraint  (Option  10  =  1).  The  second  option  cards  for 
both  inside  programs  indicate  that  the  second-order  method,  called 
generalized  Newton-Raphson,  is  used  to  minimize  the  P-function 
(Option  2=1). 

The  printout  of  results  gives  initial  information  for  the  outside 
program. 

It  indicates  that  Problem  A  is  being  solved.  It  gives  initial 
information  for  Problem  A,  noting  that  SUB  is  used.  It  solves 
Problem  A  for  the  initial  starting  point  of  the  outside  program, 
requiring  16  points,  and  indicates  that  Problem  A  is  solved. 

It  indicates  that  Problem  B  is  being  solved.  It  gives  initial 
information  for  Problem  B,  noting  that  SUB  is  used.  It  solves 
Problem  B  for  the  initial  starting  point  of  the  outside  program, 
requiring  16  points,  and  indicates  that  Problem  B  is  solved. 

It  prints  out  information  for  the  outside  program,  namely  the 
starting  point  values  of  the  constraints  found  by  solving  Problem  A 
and  Problem  B,  and  the  feasible  starting  point  to  be  used  in  the 
overall  procedure. 

It  then  prints  that  Problem  A  is  being  solved.  The  initial 
information  is  not  printed  this  time,  however,  it  having  been  by¬ 
passed  in  SUB  after  the  first  time.  Again,  Problem  A  is  solved. 

The  printout  is  given  up  to  point  4  of  this  solution.  The  remain¬ 
der,  not  included  here,  is  similar  for  more  points  of  the  outside 
program. 
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User-Supplied  Information  Cards 
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Printout  of  Results 
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MAX,  TIME*  3,6000000E*02  P*  1 .  OOOOOOOE*Ol  RATIO*  1 .OOOOOOOE^Ol  EPSILON*  1 , OOOOOOOE-OS  THET*“  1 . OOOOOOOE-OS 
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NOT  INCLUDING  TME  NON-NEGAtl VI TIES 
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NOT  1NCLU01NG  THE  NON-NEGATIVITIES 
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NOT  INCLUOING  THE  NON-NEGATIVITIeS 
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NOT  INCLUOING  THE  NON-NEGATIVITIES 
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PRDPLEH 
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NONLINEAR  PROGRAMMING  ROUTINE-SUMT  VERSION  4  SUB 
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1ST  order  estimates 
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the  current  value  of  x  15 


Lagrange  multipliers 

F«  1.036479AE*01  M  1 •0747257E*0l  G*  1 .0363509E*01  R51GHA*  0 
The  current  value  of  x  is 

1.16585UE«01  1  • 165851 1 E*01  1 .165851 1E*01  1.165851 1E*01 
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RESULTS  PRINTING  OUT  ALL  POINTS  OF  OUTSIDE  PROGRAM 

The  SUMT  program  is  modified  by  changing  RESTNT,  and  the 
INSUMT  program  is  modified  by  changing  BODYS,  CONVRX,  ESTIMS, 
FEASS,  INVERX,  OPTS,  OUTPUX,  PUNCHS,  TCHECX,  and  TIMECS  to 
eliminate  printing  of  inside  programs. 

The  user-supplied  information  cards  are  identical  to  those  of 
the  previous  section,  and  are  not  repeated  here.  The  printouts 
of  the  inside  programs  called  for  by  their  information  cards  are 
by-passed  by  the  changes  in  subroutines  mentioned  above. 

The  printout  of  the  solution  prints  initial  information  for 
the  outside  program,  for  Problem  A,  and  for  Problem  B.  From  then 
on  all  of  the  printout  pertains  to  the  outside  program.  At 
point  17  the  value  of  r  has  been  reduced  to  l.E-04,  and  the 
solution  has  been  approximately  attained.  The  program  cuts  off 
at  central  processor  time  equals  360  seconds. 
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Printout  of  Results 
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